Volume 115, Number 1, January-February 2010 

Journal of Research of the National Institute of Standards and Technology 

[J. Res. Natl. Inst. Stand. Technol. 115, 23-59 (2010)] 

Assessment of Uncertainties for the 

NIST 1016 mm Guarded-Hot-Plate Apparatus: 

Extended Analysis for Low-Density 

Fibrous-Glass Thermal Insulation 

Author's Note: This paper is based on NIST Technical Note 1606, Assessment of Uncertainties for the NIST 1016 mm Guarded-Hot-Plate 
Apparatus: Extended Analysis for Low-Density Fibrous-Glass Thermal Insulation, February 2009. Most of the content remains the same; the 
uncertainty analysis has been updated to reflect recent modifications in the heat flow imbalance study (p. 44). 



Volume 115 

Robert R. Zarr 

National Institute of Standards 
and Technology, 
Gaithersburg, MD 20899-8632 

robert.zarr@nist.gov 



Number 1 

An assessment of uncertainties for the 
National Institute of Standards and 
Technology (NIST) 1016 mm Guarded- 
Hot-Plate apparatus is presented. The 
uncertainties are reported in a format 
consistent with current NIST policy on the 
expression of measurement uncertainty. The 
report describes a procedure for determina- 
tion of component uncertainties for thermal 
conductivity and thermal resistance for the 
apparatus under operation in either the 
double-sided or single-sided mode of 
operation. An extensive example for com- 
putation of uncertainties for the single-sided 
mode of operation is provided for a 
low-density fibrous-glass blanket thermal 
insulation. For this material, the relative 
expanded uncertainty for thermal resistance 
increases from 1 % for a thickness of 
25.4 mm to 3 % for a thickness of 
228.6 mm. Although these uncertainties 
have been developed for a particular 
insulation material, the procedure and, to a 
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lesser extent, the results are applicable to 
other insulation materials measured at a 
mean temperature close to 297 K (23.9 °C, 
75 °F). The analysis identifies dominant 
components of uncertainty and, thus, 
potential areas for future improvement in 
the measurement process. For the NIST 
1016 mm Guarded-Hot-Plate apparatus, 
considerable improvement, especially at 
higher values of thermal resistance, may be 
realized by developing better control strate- 
gies for guarding that include better meas- 
urement techniques for the guard gap ther- 
mopile voltage and the temperature sensors. 

Key words: building technology; fibrous 
glass blanket; guarded hot plate; thermal 
conductivity; thermal insulation; thermal 
resistance; uncertainty. 

Accepted: November 3, 2009 

Available online: http://www.nist.gov/jres 



1. Introduction 

In October 1992, NIST officially adopted a new 
policy [1] for the expression of measurement uncertain- 
ty consistent with international practices. The NIST 
policy is based on recommendations by the Comite 
International des Poids et Mesures (CIPM) given in the 
Guide to the Expression of Uncertainty in Measurement 
[2] hereafter, called the GUM. 1 This report assesses the 



Certain commercial entities, equipment, or materials may be iden- 
tified in this document in order to describe an experimental proce- 
dure or concept adequately. Such identification is not intended to 
imply recommendation or endorsement by the National Institute of 
Standards and Technology, nor is it intended to imply that the enti- 
ties, materials, or equipment are necessarily the best available for the 
purpose. 



uncertainties for the NIST 1016 mm Guarded-Hot- 
Plate apparatus and expresses the uncertainties in a 
manner consistent with NIST policy. The uncertainty 
assessment presented herein elaborates on a previous 
effort [3] presented in 1997 for the production of NIST 
Standard Reference Material (SRM) 1450c and super- 
sedes the previous error analysis prepared by Rennex in 
1983 [4]. Technical details of the apparatus design and 
fabrication have been described previously [5-6] and, 
therefore, are only briefly presented here. 

The guarded-hot-plate method was standardized in 
1 945 after many years of effort and designated ASTM 
Test Method C 177 [7]. Essentially, the method estab- 
lishes steady-state heat flow through flat homogeneous 
slabs — the surfaces of which are in contact with ad- 
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joining parallel boundaries (i.e., plates) maintained at 
constant temperatures. The method is considered an 
absolute measurement procedure because the resulting 
heat transmission coefficients are directly determined. 
That is, the test results are not determined by ratio of 
quantities. In principle, the method can be used over a 
range of temperatures but, in this report, the mean tem- 
perature is limited primarily to 297 K (23.9 °C, 75 °F). 
This report discusses the measurement principle and 
presents a procedure for the assessment of uncertainties 
for a particular lot of low-density fibrous-glass thermal 
insulation maintained by the NIST Building and Fire 
Research Laboratory (BFRL). 

2. Reference Material 

The reference material of interest in this report is a 
low-density fibrous-glass blanket having a nominal 
bulk density of 9.6 kg • m~ 3 (0.6 lb ■ ft" 3 ). The material 
lots were manufactured in July 1980 in the form of 
large sheets (1.2 m by 2.4 m) at nominal thicknesses of 
28 mm and 81 mm. After receipt and preparation of the 
material, the National Bureau of Standards 2 announced 
in December 1980 a program [8] to provide thick 
"calibration transfer specimens" (CTS) on request 
for use in conjunction with the "representative thick- 



ness" provision of the U.S. Federal Trade Commission 
(FTC) rules published in 1979 [9] hereafter, called 
the "R-value Rule." The specimens were 610 mm 
square and were originally issued at thicknesses of 
25 mm, 75 mm, or 1 50 mm (two 75 mm specimens 
stacked). Recently, however, in order to satisfy more 
stringent energy efficiency requirements mandated in 
U.S. building codes, insulation manufacturers have 
begun requesting CTS at thicknesses up to 225 mm 
(three 75 mm specimens stacked). In accordance with 
test guidelines in the R-value Rule, measurements for 
customers are usually conducted at a mean temperature 
of 297 K and a temperature difference of either 22.2 K 
or 27.8 K (40 °F or 50 °F, respectively) across the 
specimen [9]. 

3. Steady-State Thermal Transmission 
Properties 

ASTM Practice C 1045 [10] provides a uniform 
calculation procedure for thermal transmission proper- 
ties of materials based on measurements from steady- 
state one dimensional methods such as ASTM Test 
Method C 177. Table 1 summarizes the generalized 
one-dimensional equations for thermal resistance (R), 
conductance (C), resistivity (r), and conductivity (A). 



Table 1. Steady-State One-Dimensional Thermal Transmission Property Equations 



Thermal Resistance 



Thermal Conductance 



Thermal Resistivity 
r, m-K-W -1 



Thermal Conductivity 



Equation 
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AAT 



AAT 



QL 

AAT 



Relationships 



C X 
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In 1901, Congress established the National Bureau of Standards 
(NBS) to support industry, commerce, scientific institutions, and all 
branches of government. In 1988, as part of the Omnibus Trade and 
Competiveness Act, the name was changed to the National Institute 
of Standards and Technology (NIST) to reflect the agency's broader 
mission. For historical accuracy, this report uses, where appropriate, 
NBS for events prior to 1988. 



Here, Q is the time-rate of one-dimensional heat 
flow (in units of watts, W) through the meter area of the 
guarded-hot-plate apparatus, A is the meter area of the 
apparatus normal to the heat flow direction (in units of 
square meters, m 2 ), AT" is the temperature difference 
across the specimen (in units of kelvins, K), andZ is the 
specimen thickness (in units of meters, m). As a rule, 
NIST provides value assignments and uncertainty for 
only R and, to a lesser extent, X for thermal insulation 
reference materials. Consequently, this paper presents 
uncertainty assessments only for thermal resistance (R) 
and thermal conductivity (A). 



24 



Volume 115, Number 1, January-February 2010 

Journal of Research of the National Institute of Standards and Technology 



4. Measurement Principle 

A guarded-hot-plate apparatus having appropriate 
plate temperature controllers can be operated in either a 
double sided mode or in a single-sided mode (also 
known as two-sided or one-sided mode, respectively). 
In principle, both modes of operation are covered in 
Test Method C 177; however, additional information 
on the single-sided mode is available in ASTM Practice 
C 1044 [11]. For completeness, this report presents 
both modes of operation but only the single-sided mode 
is examined in the uncertainty analysis. 

Double-Sided Mode 

Figure 1 shows the essential features of a guarded- 
hot-plate apparatus designed for operation near ambient 
temperature conditions. The plates are shown in a hor- 
izontal configuration with heat flow (Q) in the vertical 
(up/down) direction through the specimens. The appa- 
ratus is cylindrically symmetric about the axis indicat- 
ed in Fig.l. In the traditional double-sided mode of 
operation, specimens of the same material having near- 
ly the same density, size, and thickness are placed on 
each surface of the guarded hot plate and clamped 
securely by the cold plates. Ideally, the guarded hot 
plate and the cold plates provide constant-temperature 
boundary conditions to the specimen surfaces. Ideally, 
lateral heat flows (g gap and £>ed g e) are reduced to 
negligible proportions with proper guarding and, under 
steady-state conditions, the apparatus provides one- 
dimensional heat flow (Q) normal to the meter area of 



Cold plate 1 



the specimen pair. Typically, a secondary guard is 
provided by an enclosed chamber that conditions the 
ambient gas (usually air) surrounding the plates to a 
temperature near to the mean specimen temperature 
(i.e., average surface temperatures of the hot and cold 
plates in contact with the specimens). 

Under steady-state conditions, the operational defini- 
tion [10] for the mean (apparent) thermal conductivity 3 
of the specimen pair (A exp ) is 



Q 



A[(AT/L\+(AT/L) 2 ] 



(1) 



where: 

Q = the time rate of one-dimensional heat flow 
through the meter area of both specimens and, 
under ideal conditions, is equal to Q m , the elec- 
trical power input to the meter plate (W); 

A = the meter area normal to the specimen heat flow 
(m 2 ) (see Appendix A for derivation); and, 

(AT/L) { = the ratio of the surface-to-surface tempera- 
ture difference (T h - T c ) to the thickness (Z) for 
Specimen 1. A similar expression is used for 
Specimen 2. 



The thermal transmission properties of heat insulators determined 
from standard test methods typically include several mechanisms of 
heat transfer, including conduction, radiation, and possibly convec- 
tion. For that reason, some experimentalists will include the adjective 
"apparent" when describing thermal conductivity of thermal insula- 
tion. However, for brevity, the term thermal conductivity is used in 
this report. 



Guarded hot plate 
(Meter plate) 




Cold plate 2- 



Ambient gas 



1. Principle: 



T c <T h ; 



T c1 - T c2 " T c 



T c <T h ; 



T c2 ~ T c 



2. Practice: i c <- i n , i c -| 

Fig. 1. Guarded-hot-plate schematic, double-sided mode of operation — vertical heat flow. 
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For experimental situations where the temperature 
differences and the specimen thicknesses are nearly the 
same, respectively, Eq. (1) reduces to 



QL 



2AAT, 



(2) 



Using the relationship from Table 1, Eq. (2) can be 
rewritten to determine the thermal resistance of the 
specimen pair 



operation. In the single-sided mode of operation, auxil- 
iary thermal insulation is placed between the hot plate 
and the auxiliary cold plate, replacing one of the 
specimens shown in Fig. 1 . 

The auxiliary cold plate and the hot plate are main- 
tained at essentially the same temperature. The heat 
flow (Q f ) through the auxiliary insulation is calculated 
as follows [11]: 



Q'=CA{T' h -T^)~Q 



(4) 



R- 2AAT — 

Q 



(3) 



In the double-sided mode of operation, the thermal 
transmission properties correspond to an average 
temperature 7" given by T = {T h + T C )I2 . 



where the prime (') notation denotes a quantity associ- 
ated with the auxiliary thermal insulation and C" is the 
thermal conductance of the auxiliary insulation. The 
specimen heat flow (Q) is computed in the following 
equation: 



Q = Q,„-Q' 



(5) 



Single-Sided Mode 

Figure 2 shows the essential features of a guarded- 
hot-plate apparatus designed for operation near ambient 
temperature conditions in the single-sided mode of 



where Q„, is the power input to the meter plate. Values 
of Q' are typically less than 1 % of Q,„. For similar 
materials, Q from Eq. (5) is approximately one-half the 
value obtained for Q in Eq. (3) for the double-sided 
mode. 
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Fig. 2. Guarded-hot-plate schematic, single-sided mode of operation — heat flow up. 
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5. Apparatus 

Figure 3 shows an illustration of the NIST 1016 mm 
Guarded-Hot-Plate apparatus. The apparatus plates are 
typically configured in a horizontal orientation and are 
enclosed by an insulated environmental chamber that 
can be rotated ± 180°. The plates are made from 
aluminum alloy 6061 -T6. The plate surfaces in contact 
with the specimens are flat to within 0.05 mm and are 
anodized black to have a total emittance of 0.89. The 
hot plate is rigidly mounted on four bearing rods. Each 



cold plate can translate in the vertical direction for 
specimen installation and is supported at its geometric 
center by means of a swivel ball joint that allows the 
plate to tilt and conform to a nonparallel rigid sample. 
The clamping force is transmitted axially by extension 
rods that are driven by a stepper motor and a worm- 
drive gear. A load cell measures the axial force that the 
plate exerts on the specimen. The cold plates are con- 
strained in the radial direction by steel cables attached 
to four spring-loaded bearings that slide on the bearing 
rods. 




Fig. 3. NIST 1016 mm Guarded-Hot-Plate Apparatus. 
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Guarded Hot Plate 

The 1016 mm guarded hot plate is nominally 16.1 mm 
thick and consists of a meter plate 4 405.6 mm in dia- 
meter and a co-planar, concentric guard plate with an 
inner diameter of 407.2 mm. The circular gap (also 
known as "guard gap") that separates the meter plate and 
guard plate is 0.89 mm wide at the plate surface. The 
cross-sectional profile of the gap is diamond shaped in 
order to minimize lateral heat flow across the gap. The 
meter plate is supported within the guard plate by three 
stainless steel pins, equally spaced around the circum- 
ference of the meter plate, that are used to adjust the gap 
to a uniform width and maintain the meter plate in plane 
with the guard plate. Across its diameter, the meter plate 
is flat to within 0.025 mm. 

The hot-plate heater design, described previously in 
detail by Hahn et al. [12], utilizes circular line-heat 
sources located at prescribed radii. The circular line- 
heat-source for the meter plate is located at a radius of 
V2/2 times the meter-plate radius which yields a dia- 
meter of 287 mm. This location for the heater results in 
a temperature profile such that the temperature at the gap 
is equal to the average meter-plate temperature [12]. The 
heating element is a thin nichrome ribbon filament 
network, 0. 1 mm thick and 4 mm wide, electrically 
insulated with polyimide, having an electrical resistance 
at room temperature of approximately 56 Q. 



There are two circular line-heat-sources in the guard 
plate located at diameters of 524.7 mm and 802.2 mm. 
The heating elements are in metal-sheathed units, 
1.59 mm in diameter, and were pressed in circular 
grooves cut in the surfaces of the guard plate. The 
grooves were subsequently filled with a high-tempera- 
ture epoxy. The electrical resistances at room tempera- 
ture for the inner and outer guard heaters are approxi- 
mately 72 Q and 108 Q, respectively. 

Meter-Plate Electrical Power 

Figure 4 shows the electrical circuit schematic for the 
meter-plate power measurement which consists of a 
four-terminal standard resistor, nominally 0.1 Q, in 
series with the meter-plate heater. A direct-current power 
supply (40 V) provides current (/) to the circuit which is 
determined by measuring the voltage drop (V s ) across 
the standard resistor (Fig. 4) placed in an oil bath at 
25.00 °C. The voltage across the meter-plate heater (V m ) 
is measured with voltage taps welded to the heater leads 
in the center of the gap (described above). The meter 
plate power (Q m ) is the product of V m and i. 



Terminology for the 1016 mm guarded hot plate reflects current 
usage in ASTM Practice C 1043. 
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Fig. 4. Electrical schematic for meter-plate power measurement. 
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Cold Plates 

The cold plates are fabricated from 6061 -T6 
aluminum and contain channels that circulate an 
ethylene glycol/distilled water solution. Each plate is 
25.4 mm thick and consists of a 6.35 mm thick cover plate 
bonded with epoxy to a 1 9.05 mm thick base plate. The 
base plate has milled grooves 9.5 mm deep and 19.1 mm 
wide arranged in a double-spiral configuration. This 
arrangement forms a counter-flow heat exchanger, that 
is, the supply coolant flows next to the return coolant 
providing a uniform temperature distribution over the 
cold-plate surface. The temperature of each cold 
plate is maintained by circulating liquid coolant from a 
dedicated refrigerated bath regulated to within ± 0.05 K 
over a temperature range of -20 °C to 60 °C. The outer 
surfaces and edges of the cold plates are insulated with 
102 mm of extruded polystyrene foam. 

Environmental Chamber 

The environmental chamber is a large rectangular 
compartment having inside dimensions of 1 .40 m square 
by 1.60 m high supported by a horizontal axle on rota- 
tional rollers that allow the apparatus to pivot by ± 180° 
(Fig. 3). Access to the plates and specimens is permitted 
by front-and-back double-doors. Air is circulated by a 
small fan in the chamber and is conditioned by a small 
cooling coil/reheat system located within the chamber. 
N 



w 




The air temperature ranges from about 5 °C to 60 °C and 
is maintained to within ± 0.5 K by using the average of 
five Type T thermocouples located in the chamber. 

Primary Temperature Sensors 

The primary temperature sensors are small capsule 
platinum resistance thermometers (PRTs). The sensor 
construction is a strain-free platinum element supported 
in a gold-plated copper cylinder 3.18 mm in diameter by 
9.7 mm long backfilled with helium gas and hermetical- 
ly sealed. The sensors are designed for temperatures 
from 13 K to 533 K (- 260 °C to 260 °C) and the nomi- 
nal resistance is 100 Q at °C. The electrical resistance 
of each 4-wire PRT is measured with a digital multi- 
meter (DMM) that is part of an automated data acquisi- 
tion system. 

Figure 5 shows the locations of the PRTs in the cold 
plates and hot plate. The cold plate PRT is inserted in a 
3.26 mm diameter hole, 457 mm long, bored into the 
side of the cold plate (Fig. 5a). The hot plate PRT is 
located in the guard gap at an angle of 69° from the loca- 
tion where meter plate heater wires cross the guard gap 
(Fig. 5b) based on the theoretical temperature distribu- 
tion T(r, 0) determined by Hahn et al. [12] for a similar 
apparatus. The sensor is fastened with a small bracket 
on the meter side of the gap at the mid-plane of the 
plate (z = 0) as illustrated in Fig. 5c. The radius to the 
center of the PRT was computed to be 199.3 mm. 

N 



W 




Fig. Sa. Location of cold plate 
PRT (top view). 



Fig. Sb. Location of hot plate PRT in guard gap 
on meter side of guard gap (top view). 
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Meter plate 



PRT 




r 



Bracket 



Fig. 5c. Cross-section view of PRT in guard gap (guard plate not shown). 

29 



Volume 115, Number 1, January-February 2010 

Journal of Research of the National Institute of Standards and Technology 



Temperature Sensors in the Guard Gap 

The temperature difference across the guard gap 
(AT gap ) is estimated using an eight junction (4 pairs) 
Type E 5 thermopile. The thermopile was constructed 
from No. 30 AWG (American wire gauge) insulated 
thermocouple wire 0.25 mm in diameter welded in an 
argon atmosphere to form small bead junctions. The 
wire lengths were taken from spools of wire that were 
scanned using a large temperature gradient (i.e., a bath 
of liquid nitrogen) to isolate inhomogeneities in the 
wire. The wire passed from ambient to liquid nitrogen 
temperature and back to ambient; sections that gave 
thermoelectric voltages larger than 3 uV for EP wire 
and 1.7 uV for EN wire were discarded. 

Figure 6 shows the angular locations for the individ- 
ual junctions in the guard gap. The reference angle of 
0° is the location where the meter-plate heater leads 
cross the gap (the same as Fig. 5b). 



Type E is a letter designation for an ANSI standard base-metal 
thermocouple. Thermoelectric elements are designated by two letters 
where the second letter, P or N, denotes the positive or negative 
thermoelement, respectively. 



The thermocouple beads are installed in brackets 
with a thermally conductive epoxy and fastened, in 
alternating sequence, to either the meter plate or the 
guard plate similarly to the method used for the meter- 
plate PRT (Fig. 5c). Like the PRT, the junctions are 
located at the mid-plane of the hot plate (that is, z = 
in the axial direction). The EN leads of the thermopile 
depart the guard gap at an angle of 185° (as shown in 
Fig. 6) and are connected to copper leads on an iso- 
thermal block mounted inside a small aluminum enclo- 
sure. The aluminum enclosure is located inside the 
environmental chamber surrounding the apparatus 
plates. 

Temperature Control 

The three heaters in the guarded hot plate are con- 
trolled by a digital proportional, integral (PI) control 
algorithm that operates by actively controlling the plate 
temperatures. In other words, the power level is not 
fixed at a specific level which could lead to temperature 
drift. Under steady-state conditions, the meter plate 
temperature is controlled to within ± 0.003 K. 
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Fig. 6. Angular locations of Type E thermopile junctions in the guard gap (not to scale). 
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6. Measurement Uncertainty Estimation 

This section summarizes relevant uncertainty 
terminology consistent with current international 
guidelines [1-2] and presents a procedure for the esti- 
mation of measurement uncertainty based on practical 
experiences by analytical chemical laboratories [13]. 
Using this procedure, an example is given for computa- 
tion of the measurement uncertainty of the low-density 
fibrous-glass thermal insulation issued by NIST as a 
CTS. 



Terminology 

The combined standard uncertainty of a measure- 
ment result, u c (y) is expressed as the positive square 
root of the combined variance u^(y): 



u c (y) = J^c t u 2 (x t ). 



(6) 



Equation (6) is commonly referred to as the "law of 
propagation of uncertainty'" or the "root-sum-of- 
squares." The sensitivity coefficients (c,) are equal to 
the partial derivative of an input quantity (df/dXj) eval- 
uated for the input quantity equal to an input estimate 
(Xj = x,). The corresponding term, w(x,), is the standard 
uncertainty associated with the input estimate x,. The 
relative combined standard uncertainty is defined as 
follows (where y # 0): 



u c, r (y) : 



u c {y) 

\y\ 



Each u(Xj) is evaluated as either a Type A or a Type 
B standard uncertainty. Type A standard uncertainties 
are evaluated by statistical means. The evaluation of 
uncertainty by means other than a statistical analysis of 
a series of observations is termed a Type B evaluation 
[1]. Type B evaluations are usually based on scientific 
judgment and may include measurement data from 
another experiment, experience, a calibration certifi- 
cate, manufacturer specification, or other means as 
described in Refs. [1-2]. It should be emphasized 
that the designations "A" and "B" apply to the two 
methods of evaluation, not the type of error. In other 
words, the designations "A" and"B" have nothing to do 
with the traditional terms "random" or "systematic." 
Categorizing the evaluation of uncertainties as Type A 
or Type B is a matter of convenience, since both are 



based on probability distributions 6 and are combined 
equivalently. Thus, Eq. (6) can be expressed in simpli- 
fied form as: 



■\R_ 



+ u. 



(7) 



Examples of Type A and Type B evaluations are 
provided in Refs. [1-2]. A typical example of a 
Type A evaluation entails repeated observations. 
Consider an input quantity X i determined from n inde- 
pendent observations obtained under the same condi- 
tions. In this case, the input estimate x, is the sample 
mean determined from 



■ X 



1 



2X . 



The standard uncertainty, u (x,) associated with x, is 
the estimated standard deviation of the sample mean 
(where s is the standard deviation of n observations): 



u(x i ) = s(X t ) ■ 



s]n 



(8) 



The expanded uncertainty, U, is obtained by multi- 
plying the combined standard uncertainty, u c (y), by a 
coverage factor, k when an additional level of uncer- 
tainty is required that provides an interval (similar to a 
confidence interval, for example): 



U 



-- ku c (y) = ^/X c ' " 2 U ),i + X <f " 2 ( x i )b 



(9) 



The value of A: is chosen based on the desired level of 
confidence to be associated with the interval defined by 
U and typically ranges from 2 to 3 . Under a wide vari- 
ety of circumstances, a coverage factor of k = 2 defines 
an interval having a level of confidence of about 95 % 
and k = 3 defines an interval having a level of confi- 
dence greater than 99 %. At NIST, a coverage factor of 
k = 2 is used, by convention [1]. The relative expanded 
uncertainty is defined as follows (where y 2 0): 



U r 



U_ 



Note that the probability distribution for a Type B evaluation, in 
contrast to a Type A evaluation, is assumed based on the judgment of 
the experimenter. 
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For Type A evaluations, the degrees of freedom, v, is 
equal to n - 1 for the simple case given in Eq. (8). 
For the case when u c is the sum of two or more variance 
components, an effective degrees of freedom is ob- 
tained from the Welch- Satterthwaite formula as 
described in Refs. [1-2]. For Type B evaluations, V is 
assumed to be infinity. As will be shown later in this 
report, the Type B evaluation is the dominant compo- 
nent of uncertainty. Therefore, values for v are not 
necessary and are not ultimately used in determination 
of the coverage factor, k. 

Procedure 

The EURACHEM/CITAC Guide [13] provides a 
practical guide for the estimation of measurement uncer- 
tainty based on the approach presented in the GUM [2]. 
Although developed primarily for analytical chemical 
measurements, the concepts of the URACHEM/CITAC 
Guide are applicable to other fields. The primary steps 
are summarized below: 

• Specification of the mathematical process (meas- 
urement) model — clear and unambiguous state- 
ment of the measurand, i.e., Y=f(X h X 2 , ... X N ). 

• Identification of uncertainty sources — a compre- 
hensive (although perhaps not exhaustive) list of 
relevant uncertainty sources. A cause-and-effect 
diagram is a useful means for assembling this list. 

• Quantification of the components of the uncer- 
tainty sources — a detailed evaluation of the 
component uncertainties using Type A and/or 
Type B evaluations described above (for example, 
Eq. (8)) or in the GUM. 

• Calculation of the combined standard uncertain- 
ty — propagate the component uncertainties 
using the "law of propagation uncertainty" given 
in Eq. (6). 

• Calculation of the expanded uncertainty — using 
a coverage factor of k = 2, compute an interval 
for the expanded uncertainty given in Eq. (9). 



QL (Q m -AQ)L 



where: 

Qm 

AQ 



zlgap 

ft 
ft 

L 

A 
AT 



A(AT) 

(fi„-G„ 



-Of- 



■T c ) 
Q S )L 



(10) 



AT h ~T) 



power input (W) to the meter plate heater; 

parasitic heat transfer (W) from the meter 

area (defined more specifically as Q gap , ft 

and ft); 

lateral heat flow (W) across the guard gap 

(i.e., the airspace separation between the 

meter plate and guard plate shown in 

Fig. 2); 

heat flow (W) through the meter section of 

the auxiliary insulation (Fig. 2); 

error due to edge heat transfer (W) (i.e., 

from Q edge shown in Fig. 2); 

in-situ thickness of the specimen during 

testing (m); 

meter area normal to Q(m 2 ); 

specimen temperature difference (K); 

temperature of hot plate (K); and, 

temperature of cold plate (K). 



For R, the mathematical process model is given by 



R = 



A(AT) _ A{T h -T c ) 

Q Q m -^Q 

A(T h -T c ) 

Q m -Q„-Qf-Q, ' 



(ii) 



One of the major differences between Eqs. (10) and 
(1 1) is the absence of the term for specimen thickness (L) 
in Eq. (11). With regards to sign convention for heat flow 
(ft), heat gain to the meter area is assumed to be positive 
(+) and heat loss is assumed to be negative (-). 



7. Mathematical Process Model 



8. Sources of Uncertainty 



Mathematical process models are specified for 
thermal conductivity (A) and thermal resistance (R) as 
determined using the single-sided mode of operation 
(Fig. 2). For A, the mathematical process model is 
given by 



Figure 7 shows a cause-and-effect diagram that has 
been developed for A exp from Eq. (10). The cause-and- 
effect diagram is a hierarchical structure that identifies 
the main sources (shown as arrows directly affecting 
A exp ) and secondary factors (shown as arrows 
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Fig. 7. Cause-and-effect chart for A (2 levels of contributory effects). 



affecting Q, L, A, and AT") of contributory uncertainty. 
Tertiary (and additional hierarchical) factors of contrib- 
utory uncertainty are not shown in Fig. 7. In general, 
the uncertainty sources in Fig. 7 can be grouped in one 
of three major metrology categories — dimensional 
metrology for meter area (A) and thickness (Z); thermal 
metrology for temperature (7); and, electrical metrolo- 
gy for voltage (V) and resistance (Q) measurements. 
The analysis of parasitic heat losses and/or gains (A0 
requires either additional heat-transfer analyses or 
experiments (or both). 



From Fig. 7, a comprehensive, but not exhaustive, 
list of uncertainty sources is developed as shown in 
Table 2. This particular list could be applied to other 
apparatus but is most applicable to the NIST 1016 mm 
Guarded-Hot-Plate apparatus for single-sided measure- 
ments of low-density fibrous-glass blanket thermal 
insulation. Other materials, mode of operation, appara- 
tus, etc. may require a (slightly) different listing of 
sources (see, for example, the uncertainty analysis for 
NIST SRM 1450c [3]). 



Table 2. List of Uncertainty Sources for X for the NIST 1016 mm Guarded-Hot-Plate Apparatus 



1) Meter area (A) 

a) Plate dimensions 

b) Thermal expansion effects 



2) 



Thickness (L) 

a) In-situ linear position measurement system 
i) Multiple observations 

ii) System uncertainty 

b) Dimensions of fused-quartz spacers 
i) Repeated observations 

ii) Caliper uncertainty 

c) Short-term repeatability 

d) Plate flatness 

i) Repeated observations 
ii) Coordinate measuring machine (CMM) 
uncertainty 

e) Plate deflection under axial loading of cold plate 



3) Temperature difference (AT) 
a) Measurement (T h , T c ) 



i) Digital multimeter (DMM) uncertainty 
ii) PRT regression uncertainty in fit for calibration 
data 

b) Calibration of PRTs 

c) Miscellaneous sources (not shown in Fig. 7) 
i) Contact resistance 

ii) Sampling of planar plate temperature 
iii) Axial temperature variations 

4) Heat flow (Q) 

a) DC power measurement (Q m ) 
i) Standard resistor calibration 
ii) Standard resistor drift 

iii) PRT power input 
iv) Voltage measurement 

b) Parasitic heat flows (AQ) 
i) Guard-gap (Q gap ) 

ii) Auxiliary insulation (Q r ) 
iii) Edge effects (Q E ) 



The list of contributary sources of uncertainty for R 
is the same as the list given in Table 2 except the con- 



tributory source for L would be omitted, as shown in 
Eq.(ll). 
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9. Quantification of Uncertainty 
Components 

Analysis of the standard uncertainties for meter area 
(A), thickness (L), temperature difference (AT"), and 
power (Q) are presented in this section. A useful 
approach that is followed in this report is to treat each 
uncertainty component separately and evaluate the 
uncertainty component as either a Type A or Type B 
standard uncertainty [1-2]. The example presented 
here is for specimens of low-density fibrous-glass 
thermal insulation taken from the CTS lot of reference 
material in thicknesses of 25.4 mm, 76.2 mm, 152.4 mm, 
228.6 mm. The guarded-hot-plate measurements were 
conducted at a mean temperature of 297 K and a 
temperature difference of 22.2 K. The apparatus was 
operated in the single-sided mode of operation utilizing 
a specimen of expanded polystyrene foam having a 
nominal thickness of 100 mm as the auxiliary insula- 
tion (Fig. 2). 

Meter Area (A) 

The meter area is the mathematical area through 
which the heat input to the meter plate (Q) flows 
normal to the heat-flow direction under ideal guarding 
conditions (i.e., Q gap = Q e = 0) into the specimen. It is 
important to emphasize that the meter area is not the 
same as the area of the meter plate (shown in Figs. 1 
and 2). The circular meter area was calculated from 
Eq. (12) below (see Appendix A for derivation): 



^(r r f+r?)(\+aAT mp ) 2 



(12) 



c r = dA/dr = nr„(l 
c r = dA/dr t = nr t (l 

c a = dA/da = nAT, r 

dA/d(AT,J = na(r 2 



aAT mp y 
aAT mp f 

2 i „2\ 



mp ii' 1 +r, l )(\ + aAT mp ) 



C AT 



r?){\ + aATJ 



where: 



Plate Dimensions: The design gap dimensions [5] for 
the meter plate and the guard plate diameters are 405.64 
mm (15.970 in) and 407.42 mm (16.040 in), respec- 
tively. In 1994, as part of an extensive sensor 
calibration check, the meter plate was separated and 
removed from the guard plate. Using a coordinate 
measuring machine, the roundness of the meter plate 
was checked at six locations at the periphery and the 
diameter was determined to be 405.67 mm (15.971 in). 
During re-assembly, a uniform gap width of 0.89 mm 
(0.035 in) was re-established using three pin gages 
spaced at equiangular intervals between the meter plate 
and guard plate. The uncertainty of the pin gages was 
+ 0.005 mm / - 0.000 mm. Based on these check meas- 
urements, the input values for r and r„ were determined 
to be 0.20282 m and 0.20371 m, respectively, and the 
standard uncertainty for both input values was taken to 
be 0.0254 mm (0.001 in.). 

Thermal Expansion: For a, an input value of 
23.6xlO~ 6 K" 1 was taken from handbook data for 
aluminum alloy 6061-T6. The standard uncertainty 
for the value of a was assumed to be 10 % (that is, 
2.36 x 10 * K" 1 ). For tests conducted at a mean temper- 
ature of 297 K and a specimen temperature difference 
of 22.2 K, the meter plate temperature (T h ) was main- 
tained at 308 K(35 °C, 95 °F); thus, AT mp was equal to 
+ 15 K. The standard uncertainty for AT mp was deter- 
mined to be 0.086 K (and will be discussed later in the 
section on AT" uncertainty). 



r = outer radius of meter plate (m); 
r { = inner radius of guard plate (m); 
a - coefficient of thermal expansion of aluminum 

(alloy 6061-T6)(K-'); and, 
AT = temperature difference of the meter plate from 

ambient (K) = T h - 20 °C. 

The application of Eq. (6) to Eq. (12) yields 



cW{r ) + c 2 r u 2 {r i ) + c 2 a u 2 {a)+c 2 u 2 (AT ) 



u c {A): Substituting the above input estimates into 
Eq. (12), yields a meter area (A) of 0.12989 m 2 . For 
u c (A), the input estimates (x,), sensitivity coefficients 
(c,), standard uncertainties (w(x,)), and evaluation 
method (Type A or B) are summarized in Table 3. The 
last column in Table 3 provides values for c, • u (x,) to 
assess the uncertainty contribution for each input X t . 
The combined standard uncertainty u c {A) and relative 
standard uncertainty u cr {A) were determined to be 
2.4732 x 10 5 m 2 and 0.019%, respectively. This 
estimate for u c (A) is quite small near ambient tempera- 
ture but increases as T h departs from ambient condi- 
tions. 



with 
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Table 3. Summary of Standard Uncertainty Components for Meter Area (A) 



X; 



U (Xj) 



Type 



• « (*,•) 



r o 


0.20282 m 


0.63763 m 


0.0000254 m 


B 


16.20 xlO^m 2 


r i 


0.20371 m 


0.64042 m 


0.0000254 m 


B 


16.27 xlO^m 2 


a 


23.6 x 10"* KT 1 


3.8953 m 2 • K 


2.36xl0" 6 K" 1 


B 


Q.^xlO^m 2 


AT m p 


15K 


6.13 X 10~ 6 m 2 K _1 


0.086 K 


B 


O^xlO^m 2 


Thickness (L) 













In the single-sided mode of operation, the in-situ 
thickness of the specimen (Fig. 2) is monitored during 
a test by averaging four linear position transducers 
attached to the periphery of the cold plate at approxi- 
mate 90° intervals. 7 Each device consists of a digital 
readout and a slider that translates in close proximity to 
(but not in contact with) a 580 mm precision tape scale 
bonded to a precision ground plate of a low thermal 
expansion iron-nickel (FeNi36) alloy. In operation, the 
slider is excited with a pair of oscillating voltages 
which are out-of-phase by 90°. The electrical windings 
on the scale are inductively coupled with the slider and 
the resulting output signal from the scale is resolved 
and processed by the digital readout. As the slider 
follows the axial movement of the cold plate, the corre- 
sponding output signal represents the linear distance 
between the translating cold plate and the stationary hot 
plate. 

The digital readouts are reset by placing a set of four 
fused-quartz spacers of known thickness between the 
cold plate and hot plate. Fused-quartz tubing was 
selected because of its low coefficient of thermal expan- 
sion (5.5 x 10" 7 K _1 ) and high elastic modulus 
(72 GPa). The tubes have nominal inner and outer dia- 
meters of 22 mm and 25 mm, respectively. Loose-fill 
thermal insulation was placed in the tubes to suppress 
any convective heat transfer. Because the fibrous-glass 
blanket CTS is compressible, the plate separation is 
maintained during a test by four fused-quartz spacers 
placed at the periphery of the specimen at the same angu- 
lar intervals as the four linear position transducers 
described above. Four sets of spacers having lengths of 
25.4 mm, 76.2 mm, 152.4 mm, and 228.6 mm cover the 
thickness range of interest for fibrous-glass blanket CTS. 

The combined standard uncertainty for L is given by 

u c (L) = i y M 2 (Z l ) + M 2 (Z 2 ) + w 2 (i 3 ) + M 2 (I 4 )+ M 2 (Z 5 ) 
(13) 

For a two-sided test (Fig. 1), eight linear positioning devices (four 
for each specimen) determine the in-situ thickness of the specimen 
pair. 



where the sensitivity coefficients are equal to unity 
(c L . = 1) and the contributory uncertainties, identified in 
Fig. 7, are 

u(L { ) = standard uncertainty of the in-situ linear 

position measurement (m); 
u(L 2 ) = standard uncertainty of the fused-quartz 

spacers (m); 
u (L } ) = standard uncertainty of the repeatability of 

the linear position measurement (m); 
u (L 4 ) = standard uncertainty of the plate flatness 

(m); and, 
u(L 5 ) = standard uncertainty of the cold plate 

deflection under axial loading (m). 

The contributory uncertainties u (Z,) are discussed in 
detail below. 

u(L{) — In-situ Measurement: During a test, the digital 
readouts are recorded manually and the estimate for 
x(ij) is determined from the sample mean of the four 
observations. Two contributory effects comprise w(£i) : 
1) multiple observations (Type A evaluation); 
and, 2) the measurement system uncertainty (Type B 
evaluation). Thermal expansion effects of the linear 
tape scales were neglected because the iron-nickel 
(36 %) alloy has a low coefficient of thermal expansion 
and the tests are conducted near ambient conditions 
of 297 K. Equation (8) is applied to evaluate the 
Type A standard uncertainty where s is the standard 
deviation of the four transducers (n = 4). The Type B 
evaluation is the uncertainty specification stated by 
manufacturer (k= 1) of 0.005 mm. Application of 
Eq. (7) yields 



Mc (Z 1 ) = 7«>(5.0xlO- 6 ) : 



where u A varies for a particular test. Estimates for 
u (L l ) A for a test thickness of 25.4 mm are summarized 
at the end of this section (see Table 5). 
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u(L 2 ) — Spacers: Two contributory effects comprise 
u(L 2 ): 1) multiple length observations (Type A evalua- 
tion); and, 2) caliper uncertainty (Type B evaluation). 
Thermal expansion effects were neglected because 
fused quartz has a low coefficient of thermal expansion 
(5.5 x 10~ 7 K ') and the tests were conducted near 
ambient conditions of 297 K. Deformation of the 
spacers under load was also neglected because of the 
cross-sectional area of tubing and the relatively high 
value for elastic modulus. The length of each spacer 
was measured under ambient conditions with digital 
calipers and x(L 2 ) was determined from the sample 
mean of four observations. Equation (8) is applied to 
evaluate the Type A standard uncertainty where s is the 
standard deviation of the four observations (n = 4). The 
Type B evaluation assumes a uniform distribution with 
an interval of 2a [2]; thus, u B = a./y3 where a is the 
smallest length interval of the caliper. The estimates for 
u A and u B vary for each set of spacers and for the type 
of measurement calipers, respectively. Estimates for 
u(L 2 ) AB for a test thickness of 25.4 mm are summarized 
at the end of this section (see Table 5). 



u(L 3 ) — Repeatability: The short-term repeatability of 
the linear position transducers was determined from a 
series of replicate measurements. For these measure- 
ments, the digital readouts were initially set to the 
length values of each set of fused-quartz spacers placed 
between the cold plate and hot plate. The cold plate was 
lifted from the spacers and subsequently lowered in 
contact with the spacers five times to check within-day 
variation. The procedure was repeated for four 
consecutive days to check the day-to-day variation 
(20 observations total). 

The standard uncertainty for u(L 3 ) was determined 
using the Type A evaluation given in Eq. (14) [14] 



U(L 3) : 



(14) 



where s a is the standard deviation of the daily averages 
(between-day variation), s d is the (pooled) within-day 
standard deviation, and r is number of replicates per 
day (r = 5). Table 4 summarizes replication statistics 



Table 4. Summary of Replication Statistics for Uncertainty Component u (X 3 ) 



Nominal 






Within-day 


Within-day 








L 


Day 


Replicates 


Average 


Standard Deviation 


S a 


Si 


«(£ 3 ) 


(mm) 






(m) 


(m) 


(m) 


(m) 


(m) 


25.4 


1 


5 


0.0254051 


3.96 x 10 6 










2 


5 


0.0254144 


4.28x10 6 










3 


5 


0.0254156 


3.29 xl0~ 6 










4 


5 


0.0254159 


5.20x10 6 


5.12 xio" 


4.24 X Vy 6 


6.37x10^ 


76.2 


1 


5 


0.0762217 


0.70 x 10 6 










2 


5 


0.0762325 


1.93x10 6 










3 


5 


0.0762376 


1.38 x 10 6 










4 


5 


0.0762325 


3.77 x 10 6 


6.69x10 6 


2.25 X 10"* 


6.98x10 6 


152.4 


1 


5 


0.152405 


3.68 x 10~ 6 










2 


5 


0.152410 


0.70 x lo" 










3 


5 


0.152411 


3.45 X 10" 6 










4 


5 


0.152409 


2.98 x 10~ 6 


2.48 x lo" 


2.95 X lO^ 6 


3.62 X 10 -6 


228.6 


1 


5 


0.228578 


10.79 xlO" 6 










2 


5 


0.228569 


7.64 x 10 6 










3 


5 


0.228582 


2.75 x 10~ 6 










4 


5 


0.228571 


2.63 x 10 6 


6.28 x 10 6 


6.88x10 6 


8.79 XlO" 6 
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for nominal specimen thicknesses of 25.4 mm, 
76.2 mm, 152.4 mm, and 228.6 mm. Values for within- 
day average and within-day standard deviation for the 
5 replicates are given in columns 4 and 5, respectively, 
and values for s a , s d , and u(L } ) for each nominal level 
of thickness are summarized in the last three columns 
of Table 4. Note that values of u(L } ) in Table 4 do not 
appear to be correlated with L. The degrees of freedom 
(v) for Eq. (14) were determined from the Welch- 
Satterthwaite formula [1] and the value is summarized 
at the end of this section (see Table 5). 

u(L 4 ) — Plate Flatness: Two contributory effects com- 
prise u(L 4 ): 1) multiple thickness observations (Type A 
evaluation); and, 2) coordinate measuring machine 
(CMM) uncertainty (Type B evaluation). As discussed 
above, the meter plate dimensions were checked with a 
CMM in 1994. The thickness of the plate was measured 
at 32 different locations using a CMM and the estimate 
for x(L A ) was determined from the sample mean of 
32 observations. The standard deviation (s) was 
0.0131 mm and, thus, the relative flatness over the 
meter plate is (0.013 mm)/(406.4 mm) = 0.003 %. It is 
interesting to note that the flatness specification given 
in C177-04 is 0.025 % [7]. Application of Eq. (8) to 
evaluate the Type A standard uncertainty yields: 

u(L 4 ) A = 1.3 lxl(T 5 m/^32~ =2.32 xlO -"m . 

The Type B evaluation is the uncertainty specifica- 
tion (k= 1) for the CMM of 0.0051 mm. Because the 
cold plate was fabricated with the same machine finish 
as the meter plate, the cold plate flatness is assumed 
to be nearly the same as the meter plate. In this case, 
Eq. (7) becomes: 



u A L 4 ) = \I 2 ( u2 a +u l)- 

Substituting the values for the Type A and Type B 
evaluations given above yields a standard uncertainty 
for L 4 of 0.0079 mm. The value of u (L 4 ) (0.0079 mm) 
is apparatus dependent and, thus, is fixed for all values 
of specimen thickness. 

u(L s ) — Cold Plate Deflection: The potential deflection 
of the (large) cold plate under a mechanical load is 
evaluated as a Type B uncertainty using classical stress 
and strain formulae developed for flat plates. As will be 
discussed below, this approach is an approximation. 
Recall that the clamping force on the specimen and 
auxiliary insulation is transmitted axially by extension 
rods (Fig. 3). The axial force is applied over a circular 



area at the center of each plate and is assumed to be uni- 
formly distributed through a ball joint connection 
between the plate and extension rod. In the single-sided 
mode of operation, the auxiliary insulation is a rigid 
specimen of expanded polystyrene foam which sup- 
ports the hot plate (Fig. 2). For a uniform load over a 
concentric circular area of radius r, the maximum 
deflection y max at the center of the cold plate is given by 
the following formula from Ref [15]. In this case, 
simple edge support is assumed because the test speci- 
men is compressible and the plate separation is main- 
tained by edge spacers. 



u ( L 5) = y„™ = 

3W(m 2 -l) 



16k Em 



(I2m+4)a- 



m +1 



- 4rin 



a (1 m +3) r 



m +1 



(15) 



where: 

W = 

m = 

E = 

t = 



total applied load (N); 

reciprocal of Poisson's ratio (dimensionless); 

modulus of elasticity (N • m" 2 ); 

thickness of the plate (m); and, 

radius of the plate (m). 



Based on load cell measurements, a conservative 
estimate for the net applied force (W) for the cold plate 
was assumed to be 356 N (80 lbf). The plate is 1.016 m 
in diameter and 0.0254 m thick and is fabricated 
from aluminum alloy 6061-T6. The values for m, E, 
and r were taken to be (0.33) ' = 3.0, 6.9xl0 7 kPa 
(10 x 10" lb f - in 2 ), and 0.305 m, respectively. 
Substituting into Eq. (15) yields a value of 0.031 mm 
for y max , which is the dominant component of the thick- 
ness uncertainty and is essentially fixed for each level 
of specimen thickness (for constant loading). 

In general, the uncertainty due to plate deflection 
depends on the apparatus plate design (i.e., dimensions 
and material), the rigidity of the test specimen, and 
the magnitude and application of the load applied. 
The major limitations for this assessment approach 



The cold plate is not simply supported as 
assumed in Eq. (15). The plate is actually con- 
strained by the fused-quartz spacers at four loca- 
tions around the periphery of the plate. 

The cold plate is not a solid plate. As discussed 
above, the cold plate is actually a composite con- 
struction to allow the flow of coolant internally 
within the plate. 
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u c (L): Table 5 summarizes the sources, sensitivity 
coefficients (c,), uncertainty components u (£,), and the 
evaluation method (Type A or B) for a thickness of 
25.4 mm (L 254 ). As described above, the component 
uncertainties are either test dependent (u(LJ), spacer 
dependent (u(L 2 )), process dependent (m(Z 3 )), or appa- 
ratus dependent (u(L 4 )) and u(L 5 )). The final two com- 
ponents are essentially fixed for all thicknesses. 
Consequently, only the first three rows of Table 5 are 
applicable for 25.4 mm thick specimens. Application of 
Eq. (13) yields a combined standard uncertainty for 
L 25A of 0.038 mm (u cr (L) = 0.15 %). It is interesting to 
note that C 1 77-04 requires that the specimen thickness 
be determined to within 0.5 % [7]. 



Table 6 summarizes w(A)> u c(L), an d u cA^) f° r spec- 
imen thicknesses of 25.4 mm, 76.2 mm, 152.4 mm, and 
228.6 mm. As discussed above, the dominant compo- 
nent for all levels of thickness is u (L 5 ), the uncertainty 
due to potential deflection of the cold plate. As a result, 
the variation of u c (L) is small over the range of 
thicknesses. One should note that the values given in 
Table 6 are valid only for the apparatus described here- 
in. Other guarded-hot-plate apparatus would have 
different sources and values for the thickness uncertain- 
ty components. For example, the uncertainty due to 
plate flatness could be much larger if proper attention is 
not given to the plate design and fabrication. 



Table 5. Summary of Standard Uncertainty Components for 25.4 mm Thickness (L 2 s 4) 





«(*/) 


Source 


C; 


Value of u (Lj) 


Type 


«(il) 


In-situ measurement 
multiple observations 
system uncertainty 


1 


20xl0~*m 
19 X 10" 6 m 
5.0x10 6 m 


B 

A (degrees of freedom = 3) 

B (equipment specification, k= 1) 


u(L 2 ) 


Spacers (nominal 25.4) 
repeated observations 
caliper uncertainty 


1 


1.9x 10 6 m 
1.1x10 6 m 
1.5x10 6 m 


B 

A (degrees of freedom =12) 

B (a/V3 where a = 2.54 X 10" 6 m) 


u(L 3 ) 


Short-term repeatability 


1 


6.4 X 10 6 m 


A (degrees of freedom = 6.8) 


u(L 4 ) 


Plate flatness 

repeated observations 
CMM uncertainty 


1 


7.9 X 10" 6 m 
2.3xl0~ 6 m 
5.1x10 6 m 


B 

A (degrees of freedom = 31) 

B (equipment specification, k= 1) 


u(L 5 ) 


Plate deflection under load 


1 


31x10"% 


B (calculation [15]) 



Table 6. Combined Standard Uncertainty u c (L) 





(L) 


"(it) 


u(L 2 ) 


u(L 3 ) 


u(L 4 ) 


u(L 5 ) 


u c (L) 


u cr (L) 


(mm) 


(m) 


(m) 


(m) 


(m) 


(m) 


(mm) 


(%) 


25.4 


20 x 10 6 


1.9x10 6 


6.4xl0 6 


7.9x10 6 


31 x 10 6 


0.038 


0.15 


76.2 


12 x 10 6 


2.4 xl0~ 6 


7.0x10* 


7.9 x 10~ 6 


31 x 10 6 


0.035 


0.05 


152.4 


12x 10 6 


7.7 xlO -6 


3.6x10*" 


7.9 xl0~ 6 


31 X 10 6 


0.035 


0.02 


228.6 


9.6x10 6 


9.5 x 10~ 6 


8.8 x 10" 6 


7.9 x 10 "* 


31 x 10 6 


0.035 


0.02 
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Temperature Difference (AT) 

As discussed above, the primary plate temperatures 
(Figs. 1-2) are monitored during a test by computing 
temporal averages of three small capsule platinum 
resistance thermometers (PRTs) (n = 240 observations 
taken over a steady-state interval of 4 h). The uncer- 
tainty sources u (T,) for the primary temperature sensors 
are discussed in detail below. Secondary temperature 
sensors such as thermocouples and thermistors located 
in the plates, and their corresponding uncertainties, are 
not discussed because these sensors are not input quan- 
tities in the mathematical process models given in 
Eqs. (10) and (11). 

u(Ti) — Measurement: During a typical CTS test (4 h 
in duration), the electrical resistances of the PRTs are 
recorded every minute by an automated data acquisi- 
tion system (n = 240). Two major contributory effects 
comprise u(T^): 1) regression equation coefficients 
(Type A evaluation); and, 2) the measurement system 
uncertainty (Type B evaluation). (The standard uncer- 
tainty for repeated observations of AT" (Type A evalua- 
tion) was less than 0.0002 K and was neglected in 
further analyses.) 

1) For each PRT, individual observations in ohms 
(Q) were converted to temperature using a curve 
fit to the calibration data (discussed below). The 
curve-fits were obtained using a statistical plot- 
ting package from NIST The residual standard 
deviation for the fit of each set of calibration data 
was "pooled" and the resulting standard uncer- 
tainty is 0.0052 K. The degrees of freedom from 
the regression analyses were aggregated for a 
value of 15. 

2) The Type B standard uncertainty for the resist- 
ance measurement assumes a uniform distribu- 
tion with an interval 2a [2] where a was deter- 
mined from the specification of the manu- 
facturer for the digital multimeter (DMM). 
For a = 0.039 Q at the 300 Q DMM range, 
u B = a/V3 = 0.022 Q. This standard uncertainty 
in ohms was propagated using the above curve 
fit to yield a standard uncertainty for temperature 
of 0.058 K. 

u(T 2 ) — Calibration: The PRTs were calibrated by the 
NIST Thermometry Group by comparison with a stan- 
dard platinum resistance thermometer in stirred liquid 



baths. The thermometer was inserted into a test tube 
partially filled with mineral oil which, in turn, was 
placed in the calibration bath. In 1981, the thermo- 
meters were calibrated at the water triple point, 10 °C, 
20 °C, 30 °C, 40 °C, and 50 °C [4]. In 1993, the thermo- 
meters were removed from the apparatus and re- 
calibrated over an extended temperature range at 
-40 °C, °C, 40 °C, 80 °C, and 120 °C. All tempera- 
tures in the 1993 calibration were based on the 
International Temperature Scale of 1990 (ITS-90). 
Based on the expanded uncertainty (k = 2) for the cali- 
bration bath temperatures of 0.01 K (Type B evalua- 
tion), the standard uncertainty was 0.005 K (k = 1). 
Recently, the cold plate PRTs have been removed from 
their respective plates and again submitted for calibra- 
tion by the NIST Thermometry Group. These results 
will be updated when the most recent calibration and 
analysis are completed. 

u(T 3 ) — Other Small or Negligible Contributors: 
Several small or negligible contributory effects include 
the following: 1) PRT self heating/contact resistance; 
2) sampling of temperatures in the meter area (r, 9); 
and, 3) temperature variations in the axial (z) direction 
(Fig. 5c). It is difficult to quantify the uncertainties of 
these contributors by separate experiments and, in 
some cases, the uncertainties are based on theoretical 
calculations or experimenter judgment. Hence, in all 
cases, the uncertainties are Type B evaluations. 

1) PRT self-heating/contact resistance-The PRT 
excitation current is 1 mA which, for a nominal 
100 Q PRT, dissipates about 0.0001 W. For the 
meter plate PRT, a thin layer of thermally 
conductive silicone paste has been applied 
around the sensor to improve thermal contact 
(Fig. 5c). For the cold plate PRTs, the thermal 
conductance of the metal-to-air-to-metal inter- 
face between sensor and plate is estimated to be 
0.058 W • K '. Thus, the temperature rise 
(0.0001 W/0.058 W • K ') is 0.0017 K. 

2) Sampling (planar) — Rennex [4] and Siu [16] 
empirically determined the temperature profiles 
of different NIST meter plates utilizing inde- 
pendent thermopile constructions. In each exper- 
iment, the thermopiles were placed on the plate 
surfaces and a test conducted with semi-rigid 
specimens. Based on the thermopile measure- 
ments, Rennex [4] ascribed an estimate for the 
sampling uncertainty to be 0.015 K. 
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3) Axial temperature variations — A rigorous ana- 
lytical analysis by B. A. Peavy published in 
Hahn et al. [12] shows that, for typical insula- 
tions, the differences between the temperature at 
the guard gap and the average surface tempera- 
ture of the meter plate is less than 0.05 % of the 
temperature differences between the hot and cold 
plates. For a specimen temperature difference of 
22.2 K, the standard uncertainty is 0.011 K. 

u c (T): Table 7 summarizes the sources, sensitivity 
coefficients (c,), uncertainty components u (7",), and the 
evaluation method (Type A or B) for the plate tempera- 
ture. Application of Eq. (6) to the uncertainty compo- 
nents in Table 7 with c, = 1 yields a value for u c (T) of 
0.061 K (Table 7, last row). For a AT of 22.2 K, u cr (T) 
is 0.27 %. By comparison, CI 77-04 specifies an uncer- 
tainty for temperature sensors of less than 1 %. The 
dominant component for u c {T) in Table 7 is the uncer- 
tainty specification for the DMM measurement of the 
PRT electrical resistance. 

uc(AT): Recall from Eq. (10) that the specimen tem- 
perature difference (AT) was determined from the 
following equation: 



AT = (T„-T C ). 



(16) 



The application of Eq. (6) to Eq. (16) and setting 
u T = u T = u T yields 



u c (AT) = Juj 



+ U-, 



2u T . 



Substitution of u c (T) = 0.061 K (Table 7) into Eq. (17) 
yields a value for u c (AT) of 0.086 K. For AT of 22.2 K, 
« cr (A7) is 0.39 % (and for single-sided tests conducted 
(for customers) at a AT of 27.8 K, (AT) decreases to 
(0.31 %). Note that the value for u c (AT) of 0.086 K was 
used in the uncertainty assessment for the meter area (A). 

Heat Flow (Q) 

Equation (10) defines the specimen heat flow (Q) as 
the difference between the power input to the meter 
plate (Q m ) and parasitic heat losses (Q gap , Q s , and Q'). 
Ideally, in the single-sided mode of operation (Fig. 2), 
the temperatures of the guard plate, ambient gas tem- 
perature, and auxiliary cold plate are maintained such 
that the parasitic heat losses are reduced to negligible 
proportions in comparison to Q m . Thus, Q is primarily 
determined by measuring the DC voltage and current 
provided to the meter-plate heater (Q m ). The equation 
for Q m is: 



■-iV 



R 



where i is the current (VJR 3 ) measured at the standard 
resistor, and V m is the voltage drop to the meter-plate 
heater measured across the voltage taps located at the 
midpoint of the guard gap. 



Table 7. Summary of Standard Uncertainty Components for T 



»(*,) 



Source 



Value ofu(T'j) 



Type 



u(T{) 



u(T 2 ) 



Measurement (T h , T c ) 1 

DMM incertainty 1 

regression uncertainty 1 

Calibration of PRTs 1 



0.058 K 
0.058 K 
0.0052 K 

0.005 K 



B 

B (<W3 where a = 0.039 Q) 

A (degrees of freedom =15) 

B (NIST Certificate, k=V) 



u(T 3 ) 



u c (T) 



Miscellaneous 1 

Contact resistance 1 

Sampling (planar) 1 

Axial variation in plate 1 



0.019 K 
0.0017 K 
0.015 K 
0.011 K 

0.061 K 



B 

B 

B (Reference [4]) 

B (Reference [12]) 
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The application of Eq. (6) to Eq. (18) yields 



«„(&,) = -J c v" 2 (K) 2 +4" 1 (Rf +c}u\v m f 



(19) 



with 



d(Q,„)/dV s =VJR s 



c R =d(QJ/dR s = -^V m 
c Vm = d(Q,„)/dV m =VJR s 

u(Q m ) — Power Input: The contributory sources 
u (Q m ) for the meter-plate power input are discussed in 
detail below. Three contributory effects comprise 
u(Q m ): 1) calibration of the standard resistor (Type B 
evaluation); 2) PRT self-heating (Type B evaluation); 
and 3) voltage measurements for V s and V m (Type A and 
Type B evaluations). 

1) Standard resistor calibration: The 0.1 ohm 
standard resistor is a commercial, double-walled 
manganin resistor [17] manufactured in 1913. 
The resistor has been calibrated by the NIST 
Quantum Electrical Metrology Division in an oil 
bath at 25.00 °C for several years, most recently 
in 2008. Figure 8 shows the historical control 
chart for the resistor from 1977 to 2008. Since 
January 1, 1990, the NIST calibrations have been 
based on the quantum Hall effect used as the 
U.S. representation of the ohm [18]. The most 
recent calibration assigned the resistor a value 
of 0.10006957 Q and an expanded uncertainty 
(k = 2) of 0.0000005 Q. Therefore, the standard 
uncertainty was 0.00000025 Q (k= 1). 



2) PRT power input: As shown in Fig. 5c, the 
meter-plate PRT is fastened to the side of the 
meter plate. Under normal operating conditions, 
the PRT will generate a small power input to the 
meter plate due to self-heating effects. The exci- 
tation current for the meter-plate PRT is 1 mA 
which, for the nominal 100 Q PRT, dissipates a 
power of about 0.0001 W. In the worst case for a 
228.6 mm (9.0 in) thick specimen, the power 
input to the meter-plate heater is about 0.6 W 
and the PRT self-heating effect is 0.0001 W/0.6 
W = 0.02 %. Thus, the effect of PRT self-heating 
was neglected for all specimen thicknesses. 

3) Voltage measurement: Two contributory effects 
comprise the voltage measurement: 1) multiple 
observations (Type A evaluation); and, 2) the 
DMM voltage measurement uncertainty (Type B 
evaluation). During a typical CTS test (4 h in 
duration), V s and V m are recorded every minute 
by an automated data acquisition system 
(n = 240). The Type A uncertainty evaluations 
for V s and V m are included later (in Table 1 6) as 
repeated observations for the input power 
u A (Q m ). The Type B standard uncertainty for the 
voltage measurements V s and V m (Fig. 4) 
assumes a uniform distribution with an interval 
2a [2]; thus, u B = alii; where a was determined 
from the 1-year DMM specification. The DMM 
ranges for V s and V m are 30 mV and 30 V, 
respectively, and the corresponding values 
for a 30m v an d «3ov are 15-0 (J.V and 3.05 mV at 
L = 25 A mm. Therefore, u B (V s ) and u B (V m ) are 
8.7 uV and 1.76 mV, respectively. Note that cal- 
ibration checks for the DMM are conducted 
every other year; the last check was in 2008. 



Careful inspection of Fig. 8 reveals a possible 
drift in the data and its presence could be indica- 
tive of other detrimental factors affecting the 
resistor itself. Annual calibrations are now 
planned to investigate the extent of the possible 
drift. During operation with guarded-hot-plate 
apparatus, the standard resistor is immersed in an 
oil bath maintained at 25.00 °C as shown in 
Fig. 4. Because the resistor is operated at the 
same temperature as the calibration temperature 
of 25.00 °C, temperature effects during opera- 
tion were neglected. 



Table 8 summarizes the input estimates (x,), sensitiv- 
ity coefficients (c,), and standard uncertainties (w(x) 
for a CTS specimen thickness of 25.4 mm. Only Type 
B evaluation methods are included in Table 8. As stat- 
ed above, the Type A uncertainty evaluations for V s and 
V m are included (in Table 1 6) as repeated observations 
for input power u A (Q m ). Substituting the values in 
Table 8 into Eq. (19) yields a combined standard uncer- 
tainty for Q m of 0.0016 W or about 0.03 % for an input 
to the meter-plate heater of 5.15 W. The combined stan- 
dard uncertainties at other specimen thicknesses are 
summarized later in this section. 
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100076 



0.100074 



0.100072 



0.100070 



0.100068 



0.100066 



0.100064 




Fig. 8. Control chart for 0.1 Q. standard resistor, S/N 21736 



Table 8. Summary of Standard Uncertainty Components for Power Input (Q m ) at L = 25 A mm 



X, 



»<X) 



Type 



y, 


0.03 V 


169 A 


8.7 X lO^V 


B (a/VJ where a = 15.0 uV) 


R, 


0.10006957 a 


50.93 V 2 • Q~ 2 


2.5xl0" 7 n 


B (NIST Certificate, k = 1) 


V m 


17V 


0.3 A 


1.76xl0" 3 V 


B (a /V3 where a = 3.05 mV) 



u(Q gap ), u(Q'), and u(Q^ — Parasitic Heat Flows: 
Although the parasitic heat flows are reduced during 
steady-state conditions to very small values (on the 
order of 1 mW, or less), the uncertainty associated with 
each term can be large as shown below. The sources (x,) 
for the parasitic heat flows are discussed individually 
below and, later, their respective uncertainties are 
determined collectively as part of an imbalance study. 

1) Qgap — Guard gap heat flow: The model for heat 
flow across the gap developed by Woodside and 
Wilson [19] is given in Eq. (20). 



AT gap is the temperature difference across the 
guard gap. Here, the term X is the specimen ther- 
mal conductivity. The terms q„ and c are a func- 
tion of the apparatus design, specimen thickness, 
and thermal conductivity[19]. 

Empirically, Q g is determined from the 
thermopile voltage (V gap ) of an eight junction 
(4 pairs) Type E thermopile across the guard gap. 



Q' 



-S„(SnAT„y. 



■ S BapV gap (21) 



Q =(q +cX)AT 

*Z gap v ^» ' g" 



(20) 



where; q represents the heat flow directly across 
the gap; cX is the heat flow distortion in the insu- 
lation specimen adjoining the gap (Fig. 2); and 



where S gap is the heat flow sensitivity in 
W- (0.V -1 , S is the Seebeck coefficient for Type E 
thermocouples in uV- K" 1 , and n is the number of 
junction pairs. At a meter plate temperature (T h ) 
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near 308 K (35 °C), the value of S is equal to 
61.87 uV-K ' [20]; thus, the sensitivity of the 
8-junction (4 pair) thermopile is 4 x 61.87 uV- K~' 
= 247.5 liV- K '. For a DMM resolution of 
0.1 uV, the temperature resolution of the thermo- 
pile is 0.0004 °C. Under balanced control, the 
variability of the gap thermopile voltage (V gap ), 
determined from actual test data, is typically 
1.5 uV or about 0.01 K (at a 3 x standard 
deviation level). 

2) Q' — Auxiliary insulation heat flow: Equation 
(4) predicts the heat flow (Q') through the meter 
section of the auxiliary insulation, which under 
normal one-sided operation is approximately 
zero. With the exception of C", the quantities in 
the right-side of Eq. (4), namely^, T\, and T\, 
are determined from the measurement test data. 
One method for determining the thermal con- 
ductance (C") of the auxiliary insulation (in 
W • m" 2 • K" 1 ) is an iterative technique described 
in Annex Al of Practice CI 044 [11]. After the 
value of C" is obtained, the standard uncertainty 
of Q' can be determined by propagation of uncer- 
tainty in Eq. (4). An alternate method is to deter- 
mine the value of the product of C'A (that is, the 
thermal conductance per unit temperature 
(W- K ')) from the imbalance study described 
below. In this case, the standard uncertainty is 
propagated through the mathematical model 
developed for the imbalance study. 

3) Q e — Effect of edge heat transfer: In general 
terms, the edge heat flow error is the distortion 
of one-dimensional heat flow through the speci- 
men meter area due to heat flows at the periph- 
ery of the specimen. Edge effects are controlled 
by appropriate guarding in the design of the hot 
plate, limiting the specimen thickness, control- 
ling the ambient temperature at the specimen 
edge, and, if necessary, the use of edge insula- 
tion. The empirical study by Orr [21] investigat- 
ed the effects of edge insulation and changes in 
ambient temperature on edge heat flow error. A 
similar approach to determine the sensitivity of 
this error by varying the ambient air temperature 
for different specimen thicknesses was imple- 
mented as part of the imbalance study discussed 
later. 

ASTM Practice C 1043 [5] provides a theoret- 
ical analysis of edge heat loss or gain based on 



analytical solutions derived by Peavy and 
Rennex [22] for both circular and square plate 
geometries. The purpose of the analysis is to pro- 
vide the user of Practice C 1043 with design 
guidance in determining the proper diameter of 
the guard plate (Fig. 2) for control of edge heat 
loss or gain. An abbreviated version of the analy- 
sis is given below. The error (e) due to edge heat 
transfer for either geometry is given by 

2(T -T) 
e = A + BX where X =— '-^ —-T c (22) 



where T h and T c are the hot and cold plate tem- 
peratures, respectively (Fig. 2); T a is the ambient 
temperature at the edge of the specimen (Fig. 2); 
and, T m is the mean specimen temperature given 
by (T h + T c )/2. For a circular plate geometry, 
coefficients A and B are given by: 



A = %W ln and B = %W 2ii 



where 



W. 



4 (hLYyL 



" x 2 { X I b 



f{nnbl jL) 



hL 

I x (nndlyL)A I (nxd/yL) 

nnX 



The terms I and /[ are modified Bessel functions 
of the first kind of order and 1 , respectively. 
The term b is the radius of gap center; d is the 
radius of the hot plate; L is the specimen thick- 
ness; and, h is the convective film coefficient at 
the edge of the specimen. The term (hL/ X) is the 
Biot number; y 2 = X r l X z where X r and X z are the 
radial and axial thermal conductivities of the 
specimen, respectively. The geometrical mean of 
the thermal conductivities is X = (X r -X z ) m . 

The following results, for which the author is 
indebted to D. R. Flynn, retired from NIST, are 
presented in Table 9 for the following parameters: 



^ = 40 and x = ^„-r,) = 2(5) =a45a 



T h ~T. 



22.2 



Table 9 summarizes coefficients A and B, and 
the error (e) due to edge heat transfer for 
specimen thicknesses of 25.4 mm, 76.2 mm, 
152.4 mm, and 228.6 mm. 
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Table 9. Predicted values for e due to Edge Heat Transfer (Peavy and Rennex [22]) 



L 


A 


B 




e 




(mm) 






X= -0.450 


x=o 


X= +0.450 


25.4 

















76.2 





0.0000021 











152.4 


0.0000021 


0.0022826 


-0.0010 





+ 0.0013 


228.6 


0.0002111 


0.0266257 


-0.0118 


0.0002 


+ 0.0122 



The error (e) for X = is essentially zero up 
to thicknesses of 228.6 mm. For an ambient 
temperature difference of ± 5 K from the mean 
specimen temperature (X= + 0.450), predicted 
values of s become significant at thickness of 
152.4 mm. 

tsQ — Imbalance study: A series of imbalance tests 
were conducted to investigate empirically the effects of 
moderate temperature differences on Q mp , Q', and Q e . 
An imbalance test, as the name implies, is an operating 
condition in which a parameter (in this case, tempera- 
ture difference) is purposely imbalanced (from zero) 
such that a parasitic heat flow is enhanced (i.e., the 
magnitude of the heat flow is either increased or 
decreased). In this study, imbalance tests are conducted 
for temperature differences across: l)the guard gap 
(AT gap or, in this case, V gap as shown in Eq. (21)); 
2) auxiliary insulation (T' h - T' c ); and, 3) the mean 
specimen temperature at the specimen edge and the 
ambient air temperature (T m - T a ). 

These parameters were varied following an ortho- 
gonal experimental design illustrated in Fig. 9, where 
x i = V gap ; x 2 = T' h - T' c ; and x 3 = T m - T a . The test plan 
for the imbalance tests is based on a 2 3 full factorial 



design meaning that the three factors are each varied at 
one of two levels shown in Yates order in the adjoining 
design matrix (with levels "coded" +1 for a high setting 
and -1 for a low setting). One advantage of an ortho- 
gonal design is that any interactions between factors 
(xfc 2 j x i x i j X 2 X 3 > an d X \ X 2 X }) can be detected; that is not 
possible for an experimental design that varies "one- 
factor-at-a-time." Note that the experimental design 
given in Fig. 9 also contains a center point (#9), 
that is, a balanced test point where all imbalance 
settings have been set equal to zero. 

For a 2 3 full factorial design, the useful underlying 
mathematical model would be 

y = P + W i x l +P 2 x 2 +p } x i +p l2 x 1 x 2 

+ Pu X l X i + Pl3 X 2 X 3+Pm X l X 2 X 3 )l 2 (23) 



where: 

Xj 

A 

Ay 



"coded" level (+1 or-1) for factor /; 
main effect for coded x,; 
two-factor interaction effect for coded x{, and, 
three-factor interaction effect for coded x,. 



7 


t x 2 




s 


/[ 




/ 




| 

- 

: 

/ 


*3 

/ 










6 




/ 





2 s Full Factorial Design Matrix - Yates Order 


# Xl X3 Xi 


1 - - - 


2 


3 - + - 


4 + + - 


5 - \ + j 


6 + - + 


7 — + + 


8 + + + 


9 



Fig. 9 Full-factorial experimental design for 3 factors, 2 levels. 
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Table 10 summarizes the input settings for the nine 
test conditions following the procedure given in Fig. 9. 
For tests #1 through #8, the "non-coded" settings for x,, 
x 2 , and x 3 are ± 50 (j,V, ± 0.5 K, and ± 5 K, respectively. 
The input values for the test settings were selected with 
the objective of providing adequate responses in the 
parasitic heat flows. Test #9 is a balanced condition 
where each parameter is set to zero and the correspon- 
ding parasitic heat flows consequently approach zero. 
Note that Test #9 is an independent check test for the 
original test measurement. 

Table 10. Nominal Settings for Imbalance Study in Yates Order 



gap 



n-n 

(K) 



T -T 

m a 

(K) 



1 


-50 


-0.5 


-5 


2 


150 


-0.5 


-5 


3 


-50 


+0.5 


-5 


4 


150 


+0.5 


-5 


5 


-50 


-0.5 


15 


6 


150 


-0.5 


15 


7 


-50 


+0.5 


15 


8 


150 


+0.5 


+5 


9 












The actual test sequence (not shown in Table 10) was 
randomized in order to minimize the introduction of bias 
in the results. The guarded-hot-plate imbalance tests 
were conducted with low-density fibrous-glass CTS 
specimens having thicknesses of 25.4 mm, 76.2 mm, 
152.4 mm, and 228.6 mm (n = 36 data points) at a mean 
temperature of 297 K and a AT^of 22.2 K across the spec- 
imen. The apparatus was operated in the single-sided 
mode of operation utilizing a specimen of expanded 
polystyrene foam having a nominal thickness of 100 mm 
as the auxiliary insulation (Fig. 2). 

Table 11 summarizes the 36 values for Q m , AT, 
V g ap, T\ - T' c ; and, T m - T a at specimen thickness of 25.4 
mm, 76.2 mm, 152.4 mm, and 228.6 mm. For a particu- 
lar specimen thickness, values of Q m varied considerably, 
depending on the imbalance settings for x l = V' x 2 = 



T' h -T': and, 



T„, The balanced condition 



provided in Test #9 was used to establish a baseline value 
for analyses of the data. At the balanced condition (Test 
#9), values of Q m decreased from about 5.1 W to 0.6 W 
as the specimen thickness increased from 25.4 mm to 
228.6 mm. This 8 fold decrease in Q m is important and 
will be treated later. 



Table 11a. Imbalance Test Data in Yates Order 





25.4 mm 


76.2 mm 




Q,n 


AT 


V 

gap 


n-n 


T -T 

m a 


Q m 


AT 


V 

gap 


n-K 


T -T 
m a 


it 


(W) 


(K) 


(UV) 


(K) 


(K) 


(W) 


(K) 


(HV) 


(K) 


(K) 


1 


4.9701 


22.22 


^19.96 


-0.501 


-5.00 


1.6137 


22.22 


-49.96 


-0.499 


-5.00 


2 


5.2257 


22.22 


+50.02 


-0.505 


-5.00 


1.8798 


22.22 


+49.99 


-0.495 


-5.00 


3 


5.0159 


22.22 


-49.97 


+0.499 


-5.00 


1.6621 


22.22 


-50.00 


+0.499 


-5.00 


4 


5.2768 


22.22 


+50.03 


+0.502 


-5.00 


1.9285 


22.22 


+50.00 


+0.497 


-5.00 


5 


4.9589 


22.22 


-49.94 


-0.503 


+5.00 


1.6165 


22.22 


-50.02 


-0.500 


+5.00 


6 


5.2153 


22.22 


+50.01 


-0.498 


+5.00 


1.8829 


22.22 


+49.96 


-0.500 


+5.01 


7 


5.0067 


22.22 


^19.96 


+0.506 


+5.00 


1.6644 


22.22 


-49.98 


+0.497 


+5.00 


8 


5.2652 


22.22 


+50.05 


+0.501 


+5.00 


1.9306 


22.22 


+50.03 


+0.497 


+5.00 


9 


5.1189 


22.23 


0.06 








1.7719 


22.22 


0.00 









Table lib. Imbalance Test Data in Yates Order 





152.4 mm 


228.6 mm 




Q,n 


AT 


V 

gap 


T'h-T" o 


T -T 

m a 


Q m 


AT 


V 

gap 


n-n 


T -T 

m a 


# 


(W) 


(K) 


(HV) 


(K) 


(K) 


(W) 


(K) 


(HV) 


(K) 


(K) 


1 


0.7276 


22.22 


-49.98 


-0.500 


-5.00 


0.4425 


22.22 


-49.99 


-0.502 


-5.00 


2 


0.9978 


22.22 


+49.99 


-0.502 


-5.00 


0.7225 


22.22 


+49.97 


-0.499 


-5.00 


3 


0.7752 


22.22 


^19.98 


+0.495 


-5.00 


0.4915 


22.22 


-49.91 


+0.499 


-5.00 


4 


1.0461 


22.22 


+49.98 


+0.493 


-5.00 


0.7714 


22.21 


+50.03 


+0.503 


-5.00 


5 


0.7269 


22.23 


-50.00 


-0.500 


+5.00 


0.4397 


22.22 


-49.99 


-0.503 


+5.01 


6 


0.9975 


22.22 


+50.04 


-0.506 


+5.00 


0.7200 


22.22 


+50.07 


-0.492 


+5.00 


7 


0.7757 


22.21 


-50.06 


+0.497 


+5.01 


0.4897 


22.22 


-49.97 


+0.500 


+5.00 


8 


1.0457 


22.22 


+50.02 


+0.500 


+5.00 


0.7660 


22.22 


+49.96 


+0.495 


+5.01 


9 


0.8852 


22.22 


+0.01 








0.6042 


22.21 


+0.03 
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Values of Q m from Table 11 and coded values (±1) of 
x b x 2 , and x 3 were input into Eq. (23) and the least- 
squares estimated effects (/$,-, ji ip and jS, /fc ) were calculat- 
ed using a NIST statistical graphical analysis program 
[23] that employed Yates algorithm [24]. Table 12 sum- 
marizes whether the estimated effects (/$,-, jS, 7 , and Py k ) are 
statistically significant at the 5 % level or the 1 % level. 
In all cases, factors x l and x 2 were significant. Somewhat 
surprisingly, however, factor x 3 was determined to be 
insignificant forZ = 152.4 mm and 228.6 mm. The effects 
estimates for all interactions (x { x 2 , x { x 3 , X2X3, and x l x 2 x i ) 
across all thicknesses were insignificant. Based on these 
results, the (non-coded) data in Table 11 were 
subsequently fitted to a simplified model (discussed 
below). 



Table 12. Statistical Significance for Estimated Effects for 
Imbalance Study 



data from Table 1 1 were fit to an empirical model for AQ 
given in Eq. (24) below, 





Factor 


25 A mm 


76.2 mm 


152.4 mm 


228.6 mm 


X 


** 


** 


** 


** 


x 2 


* 


** 


** 


* 


x i 


* 


* 


— 


— 


x \ x 2 


— 


— 


— 


— 


X \ X T, 


— 


— 


— 


— 


XjX 3 


— 


— 


— 


— 


x l x 2 x 3 


— 


— 


— 


— 



** Statistically significant at the 1 % level 
* Statistically significant at the 5 % level 

The parasitic heat flows from the meter area, expressed 
mathematically in Eq. (10), represent lateral heat flows 
due to temperature differences across the guard gap 
and the auxiliary insulation, and due to specimen edge 
effects. Based on the results from Table 12, the 



^ z^tft ; *~>m 



where 

Qm, 



■ay gap +a 2 {T;-T;) + a,{T n -T) 

(24) 



Q: 



K 

T 



the power input to the meter-plate heater for 

test i = 1 to 9; 

the power input to the meter-plate heater for 

test i = #9 (balanced condition); 
Uj = regression coefficients (/' = 1,2,3); 
Vgap = guard gap voltage corresponding to AT gap 

(HV); 

hot plate temperature (K); 

auxiliary plate temperature (K); 

mean specimen temperature (K) = (T h + T c ) 

12; and, 
T a = ambient gas temperature near the edge of the 

specimen (K). 

The presence of an offset coefficient a was consid- 
ered for the model given in Eq. (24) but, because the 
term is predicted to be nearly zero from theory (Eq. 22), 
the term is not included here. 

Table 13 summarizes estimates and approximate 
standard deviations determined by multiple variable 
linear regression for coefficients, a b a 2 , and a } as a 
function of specimen thickness. The "goodness of fit" 
denoted by the residual standard deviation (RSD) in the 
last column is equal to, or less than, 3 mW (Type A eval- 
uation). If an offset coefficient a is incorporated in 
Eq. (24), the RSDs for the fits are about 1 mW, or less. 
Regression coefficients a b a 2 , and a 3 are discussed 
below. 



Table 13. Estimates and Standard Deviations for a t , a 2 , and a 3 in Eq. (24) 



L 
(mm) 


01 1 
(W • uV~ ) 


■s(ai) 
(W-uV "') 


a 2 
(WK ) 


s(a 2 ) 
(W-K -1 ) 


° 3 1 
(W • K ) 


s(a 3 ) 

(W-K ') 


RSD* 
(W) 


25.4 


+0.002579 


2.15x10 5 


0.04846 


2.14x10 3 


+0.001072 


2.15 x 10 4 


0.0030 


76.2 


+0.002663 


4.59x10 6 


0.04833 


4.61 Xl0 4 


-0.000269 


4.59x10 5 


0.0006 


152.4 


+0.002705 


1.36x10 5 


0.04836 


1.36x10 3 


+0.000028 


1.36xl0 4 


0.0019 


228.6 


+0.002790 


1.24x10 5 


0.04856 


1.24x10 3 


+0.000309 


1.24 x 10 4 


0.0017 



*Residual standard deviation for fit 
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1) a x : Estimates for a, in Table 13 represent the 
heat flow sensitivities (S gap ) in W-|xV _1 for the 
4-pair thermopile as defined in Eq. (21). At 
25.4 mm, the a x estimate is consistent with 
previous results 8 obtained for SRM 1450c, 
Fibrous Glass Board [3], Note, however, that 
estimates fo a x increase 8.2 % from 25.4 mm to 
228.6 mm indicating that, for a fixed tempera- 
ture difference, the heat-flow across the guard 
gap increases with the specimen thermal resist- 
ance. From Eq. (20) and Table 1, the ratio of gap 
heat flow and specimen heat flow is 

Q„ (g +cX)AT„ 



Q 



AAT/R 



or 



e. 



Q 



1 AT 
A AT 



(q o R+cL). 



(25) 



Therefore, for a particular apparatus (i.e., fixed 
A) and a given temperature imbalance (AT gap ), 
the ratio of lateral and specimen heat flows 
increases as the specimen resistance (R) and 
thickness (L) increase if the other factors (name- 
ly AT) remain unchanged. Consequently, if one 
assumes a constant coefficient for a, based on L 
equal to 25.4 mm, an error in uncertainty is 
introduced for tests conducted at greater 
thicknesses. 



values obtained for the thermal conductance 
(C'J computed from Eq. (26) below: 



C' = C'A ■■ 



k'A 

IT 



(26) 



Substitution of £'=0.0373 W-rrr'-K ' obtained 

from independent guarded-hot-plate tests at 

311 K, A = 0.12989 m 2 at 308 K, and L' = 



0.1013 m in Eq. (26) yields a value 
0.0478 W-K ' which is within 1.3 % of a, . 



of 



3) a 3 : Although the estimates for a 3 in Table 13 
represent a thermal conductance per unit temper- 
ature, it is more useful to express these estimates 
as an edge heat flow error (e). By setting 
V„=0, (T'„-T' c ) = 0, and dividing by Q mo , 
Eq. (24) becomes 



Afi 



Q -Q 



z~ m 



-(T m -T). 



(27) 



Table 14 summarizes values of e computed from 
Eq. (27) for (T m - T a ) equal to K and ± 5 K and 
values of e reproduced from Table 9 for L rang- 
ing from 25.4 mm to 228.6 mm. Theoretical 
values of e are based on calculations from 
Eq. (22) for the same values of temperature 
imbalances. There are two general conclusions 
from the results presented in Table 14: 



2) a 2 : Estimates for a 2 in Table 13 represent the 
thermal conductance per unit temperature (C' a ) 
of the auxiliary specimen (100 mm thick poly- 
styrene foam) at 308 K (35 °C). The average 
value for Table 13 is -0.04843 W-K ' and range 
is 0.00022 W-K ' (or 0.5 %). As a check, the 
values for a x were compared with computed 



At specimen thicknesses less than or equal 
to 76.2 mm, the measured effect is small 
but greater than predicted by theoretical 
analysis. The fact that the empirical results 
are not monotonic with L suggests that the 
observed variations for 152.4 mm and 
below are due to experimental variations 
of factors other than L. 



Rennex [4] obtained a value of approximately 0.00057 W-uV~ 
for an 18-pair Type E thermopile used in the guard gap. In 1994, the 
18-pair thermopile was replaced, as part of upgrade and regular 
maintenance, with the 4-pair Type E thermopile presented in this 
report. Note that the difference in sensitivities between the 4-pair and 
18-pair thermopiles is 4.5, roughly the same ratio observed between 
the value reported by Rennex [4] and the values for a x reported in 
Table 13 (0.0026/0.00057 = 4.56). 



• At specimen thicknesses greater than or 
equal to 152.4 mm, the measured effect is 
much smaller than predicted by theoretical 
analysis. 

Further work is recommended to investigate the 
differences between edge heat flow errors com- 
puted from empirically-derived models and from 
theoretical based models. 
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Table 14. Comparison of Empirical and Theoretical Values of Edge Heat Flow Error (e) 



L 

(mm) 

25.4 

76.2 

152.4 

228.6 



Empirical e (Equation 27)) 

-T a = -5K T m -T a = 0K r m -T fl = +5K 

-0.0010 +0.0010 

+ 0.0008 - 0.0008 

-0.0001 +0.0001 

- 0.0025 + 0.0025 



Theoretical e (Table 9) 
X=- 0.450 X=0 X= + 0.450 






-0.0010 
-0.0118 


















+ 0.0013 


0.0002 


+ 0.0122 



u c (AQ) — Standard uncertainty for parasitic heat flows: 



Substitution of x. 



T\ - T' 



and, 



1 ' gap J A 2 i h 

T m - T a for the quantities in Eq. (24) and applica- 



tion of Eq. (6) yields 
u e (AQ) 



2 2 , x , 2 2 , x . 2 2/ x 

<c a u (a i ) + c x u (Xjj + c^u (a 2 ) 



+c x u (x 2 ) + c a u (a 3 ) + c x3 M (x 3 ) 



(28) 



with 

c a =d(AQ)/daj = x . 
c x = d(AQ) / dxj = a j . 

Table 15 summarizes input values for Eq. (28) and 
the corresponding u c (AQ) for specimen thicknesses of 
25.4 mm, 76.2 mm, 152.4 mm, and 228.6 mm. The 
input values for u (a) and a, were obtained from 
Table 13. Under steady-state test conditions, the input 
estimates for x ; - are nearly zero (Table 1 5). The standard 
uncertainties for x, (AT gap ), x 2 (T' h — T' c ), and 
x 3 (T m - T a ) were estimated to be ± 0.01 K, which corre- 
sponds to ± 2.48 uV for V gap , ± 0.086 K and ± 0.5 K, 
respectively. Note that the resulting values of u c (AQ) 
are nearly constant across all levels of specimen 
thickness. 

u c (Q) — Standard uncertainty for specimen heat flow: 
Recalling the heat balance for the meter area given in 



Eqs. (10) and (11), the specimen heat flow (Q) was 

determined from the following equation: 

(29) 

Application of Eqs. (6) and (7) to Eq. (29) yields 



u c Q = ju 2 A (Q m ) + ul{Q m ) + ul{AQ) (30) 

where the sensitivity coefficients are equal to unity and 
the contributory uncertainties are 

u a (Qm) = standard uncertainty of repeated power meas- 
urements (n = 240) during a test (W); 

u b (Qm) ~ standard uncertainty of the meter-plate 
power input — see Table 8 (W); and, 

u c (AQ) = combined standard uncertainty for parasitic 
heat flows computed in Eq. (28) — see Table 
15 (W). 

Table 1 6 summarizes the estimates for the combined 
standard uncertainty u c (Q) across all thicknesses. The 
dominant source by, in many cases, two orders of mag- 
nitude is the uncertainty in the parasitic heat flows 
which, in turn, is due primarily to the uncertainty in the 
gap thermopile voltage and temperature difference 
across the auxiliary insulation. The relative combined 
standard uncertainty u cr (Q) increases considerably 
with specimen thickness from 0.17% at 25.4 mm to 
1.4% at 228.6 mm. 



Table 15. Estimates for u c (AQ) 



L X\ 


"(fll) 


a, 


"(*i) 


x 2 


u(a 2 ) 


a 2 


u(x 2 ) 


x 3 


u(a } ) 


ai , 


u(x 3 ) 


uMQ) 


(mm) (uV) 


(W • itV' 1 ) 


(W • uV" 1 ) 


(K) 


(K) 


(W • K" 1 ) 


(W • K"') 


(K) 


(K) 


(W • K" 1 ) 


(W-K -1 ) 


(K) 


(W) 


25.4 0.01 


2.15 x 10 5 


+0.002579 


2.48 


0.005 


2.14x10 3 


0.04846 


0.086 


0.004 


2.15 X 10 4 


+0.001072 


0.5 


0.0087 


76.2 0.02 


4.59 x 10 6 


+0.002663 


2.48 


0.003 


4.61 x 10 4 


0.04833 


0.086 


0.003 


4.59x10 5 


-0.000269 


0.5 


0.0082 


152.4 0.01 


1.36 x 10 5 


+0.002705 


2.48 


0.002 


1.36x10 3 


0.04836 


0.086 


0.004 


1.36 Xl0~ 4 


+0.000028 


0.5 


0.0088 


228.6 0.01 


1.24 x 10 5 


+0.002790 


2.48 


0.002 


1.24 xlO 3 


0.04856 


0.086 


0.002 


1.24x10^ 


+0.000309 


0.5 


0.0083 
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Table 16. Combined Standard Uncertainty u cl Q) 





L 
(mm) 


Q 

(W) 


«a(QJ 

(W) 


U B (Q,n) 

(W) 


« C (A0 
(W) 


«„(G) 

(W) 


M6) 

(%) 


25.4 


5.1452 


0.0006 


0.0016 


0.0087 


0.0089 


0.17 


76.2 


1.8032 


0.0004 


0.0007 


0.0082 


0.0082 


0.46 


152.4 


0.8761 


0.0003 


0.0004 


0.0088 


0.0088 


1.0 


228.6 


0.6112 


0.0003 


0.0003 


0.0083 


0.0084 


1.4 



10. Calculation of Combined Uncertainty 



The standard uncertainties for meter area (^4). thick- 
ness (L), temperature difference (AT), and power (Q) 
calculated in the previous section are combined using 
the law of propagation uncertainty given in the GUM 
(Eq. (6)). The combined standard uncertainties (u c , u cr ) 
for X exp , and R are presented for specimens of low- 
density fibrous-glass thermal insulation at thick- 
nesses of 25.4 mm, 76.2 mm, 152.4 mm, and 
228.6 mm. The expanded uncertainties (U, U r ) for X exp 
and R are also presented for the same thicknesses using 
a coverage factor of k = 2. 

Combined Standard Uncertainty and Expanded 
Uncertainty for X ap 

Application of Eq. (6) to Eq. (10) yields the com- 
bined standard uncertainty for X exp : 



with 



u c (A exp ) = JV Q u 2 ( + c\ u c (A) + «4 u c ( A T) 



(31) 



C Q 


U'^xp' "^ 


A(AT) 


c L = 


3A exp /3Z 


Q 

A(AT) 


C A = 


dA exp /dA 


-QL 

A 2 (AT) 


C AT = 


dPi exp /dAT 


-QL 

A(AT\ 2 



Table 1 7 summarizes the input estimates (x,), sensitivi- 
ty coefficients (c,), standard uncertainties (u c (x,)) for 
Q, L, A, and AT" at thicknesses of 25.4 mm, 76.2 mm, 
152.4 mm, and 228.6 mm. The last two columns of 
each table provide the absolute and relative contribu- 
tions of each component uncertainty. The last two rows 
of each table provide A exp , u c (X exv ), u cr (X exv ), U, and U r 
using a coverage factor of k = 2. The estimate for A exp 
and corresponding bulk density (p) are also given on 
the last row of each table. 



Table 17a. Combined Standard Uncertainty for Thermal Conductivity (A exp ) for L = 25.4 mm 



Q 

L 

A 

AT 

y = 4x P 



5.1452 W 
0.0254 m 
0.12989 m 2 
22.22 K 
= 0.0454 W • m" 1 ■ K"'(p = 9.4 kg • m" 3 ) 



0.0089 W 
3.8 X 10 5 1 
2.47x10"-' 
0.086 K 



u c(Kxp) 

U(k = 2) 



\c, u\ 


(c, u)ly 


(W • m • K ) 


(%) 


0.00008 


0.17 


0.00007 


0.15 


0.00001 


0.02 


0.00017 


0.39 


0.00020 


0.45 


0.00041 


0.9 
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Table 17b. Combined Standard Uncertainty for Thermal Conductivity (A e ) for L = 76.2 mm 



















\c, u\ 


(c; u)ly 


*r 


*t 




c i 


u(x,) 






(W-m-K" 1 ) 


(%) 


Q 


1.8032 W 




0.02641 nf'-K" 1 


0.0082 W 






0.00022 


0.46 


L 


0.0762 m 




0.62466 W-m" 2 -K"' 


3.5xl0" 5 m 






0.00002 


0.05 


A 


0.12989 m 2 




-0.36657 W-m" 3 -K"' 


2.47xl0" 5 m 2 






0.00001 


0.02 


AT 


22.22 K 




-0.00214 W- nf'-K" 2 


0.086 K 






0.00018 


0.39 


y = Kx P = 


= 0.0477 W-m" 1 K" 1 (p = 


= 9.0 kg 


m' 3 ) 




u c (X, 
U(k = 


2) 


0.00029 
0.00057 


0.60 

1.2 



Table 17c. Combined Standard Uncertainty for Thermal Conductivity (A ) fori = 152.4 mm 

















\C, U\ 


(c t u)ly 


*i 


*; 






c i 


u(x,) 




(W ■ m" 1 • K" 1 ) 


(%) 


Q 


0.8761 W 






0.05280 m" 1 • K" 1 


0.0088 W 




0.00046 


1.0 


L 


0.1524 m 






0.30356 W-m" 2 -K"' 


3.5xl0" 5 m 




0.00001 


0.02 


A 


0.12989 m* 


1 




-0.35614W-m" 3 -K"' 


2.47xl0" 5 m 2 




0.00001 


0.02 


AT 


22.22 K 






-0.00208 W-m"'-K" 2 


0.086 K 




0.00018 


0.39 


y = Kx P = 


= 0.0463 W- nf'-K" 


1 (P- 


= 9.6 kg- 


m" 3 ) 




u c(Kxp) 


0.00049 


1.1 



U(k = 2) 0.00099 



2.1 



Table 17d. Combined Standard Uncertainty for Thermal Conductivity (A exp ) for L = 228.6 mm 



















\c, u\ 


(c,- u)ly 


x, 


*< 




c i 


"(*,) 






(W ■ m • K ) 


(%) 


Q 


0.6112 W 




0.07920 m^-K" 1 


0.0084 W 






0.00066 


1.4 


L 


0.12286 m 




0.21178 W-m" 2 -K _1 


3.5xl0" 5 m 






0.00001 


0.02 


A 


0.12989 m 2 




-0.37270 W-m" 3 -K"' 


2.47xl0" 5 m 2 






0.00001 


0.02 


AT 


22.22 K 




-0.00218 W-m" 1 K" 2 


0.086 K 






0.00019 


0.39 


y = 4x P = 


= 0.0484 W-m" 1 -K"'(p = 


= 8.7kg- 


m" 3 ) 




U(k = 


:xp> 

2) 


0.00068 
0.00137 


1.4 

2.8 



Using the values summarized in Table 17, Fig. 10 
plots the relative combined standard uncertainty for A 
and individual uncertainty components for Q, L, A, and 
AT" as a function of L. The graphical analysis clearly 
shows that the uncertainty contribution from Q is the 
major component of uncertainty, especially at thick- 
nesses of 75 mm and above. At 25 mm, the major com- 
ponent is the standard uncertainty of AT" which is fixed 
at 0.39 % across all levels of L. The uncertainty contri- 
bution due to L is largest at 25 mm and declines dra- 
matically, as expected, as L increases. The uncertainty 
contribution due to A is nearly zero for all thicknesses. 



Figure 11 plots the estimates for thermal conduct- 
ivity (A) of the fibrous-glass blanket CTS given in 
Table 17 as a function of specimen thickness (L). The 
error bars shown in Fig. 1 1 are equal to the expanded 
uncertainties U(k = 2) given in Table 17. The data in 
Fig. 11 show a "thickness effect" which was the basis 
for the R-value Rule [9], and the subsequent develop- 
ment of the NIST CTS lot of material. The data in 
Fig. 1 1 also reveal a material variability factor for the 
NIST CTS lot of material. 
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Fig. 10. Combined standard uncertainty and individual components for X. 
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Fig. 11. Thermal conductivity measurements of Fibrous-Glass Blanket NIST CTS as a function of thickness 



Combined Standard Uncertainty and Expanded 
Uncertainty for R 

Application of Eq. (6) to Eq. (11) yields the com- 
bined standard uncertainty for if: 



250 



u c( R ) = V4"c (0 + C X 2 ( A ) + c\ T u]{AT) 



(32) 



with 




c g = dR/dQ = 


A 2 (AT) 

Q 2 


c A = dR/dA = 


AT 
~Q 


c AT = dR/dAT = 


A 



Q 
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Table 18 summarizes the input estimates (x,), sensi- 
tivity coefficients (c,), standard uncertainties (m c (x,)) 
for Q, A, and AT" at thicknesses of 25.4 mm, 76.2 mm, 
152.4 mm, and 228.6 mm. The last two columns of 



each table provide the absolute and relative contribu- 
tions of each component uncertainty. The last two rows 
of each table provide R, u c (R), u cr (R), U and U r using 
a coverage factor of k = 2. 



Table 18a. Combined Standard Uncertainty for Thermal Resistance (R) for L = 25.4 mm 















\c, u\ 


(c; u)ly 


x, 


x i 


c i 


"(*,) 




(m 2 -K-W _1 ) 


(%) 


Q 


5.1452 W 


-0.10901 m 2 KW 2 


0.0089 W 




0.00097 


0.17 


A 


0.12989 m 2 


4.31811 K-W -1 


2.47xl0" 5 m 2 




0.00011 


0.02 


AT 


22.22 K 


0.02524 m 2 -W ' 


0.086 K 




0.00218 


0.39 


y = R- 


= 0.560 m 2 K-W ' 






u c (R) 

U(k = 2) 


0.0024 
0.0048 


0.43 
0.9 



Table 18b. Combined Standard Uncertainty for Thermal Resistance (R) for L = 76.2 mm 













|C; u\ 




(c,- u)ly 


*i 


x i 


c i 


»(■*;) 




(m 2 -K-W 


') 


(%) 


Q 


1.8032 W 


-0.88782 m 2 • K • W 2 


0.0082 W 




0.00732 




0.46 


A 


0.12989 m 2 


12.325 K-W" 1 


2.47xl0" 5 m 2 




0.00030 




0.02 


AT 


22.22 K 


0.07204 m 2 -W _1 


0.086 K 




0.00621 




0.39 


y = R-- 


= 1.60 m 2 KW ' 






u c (R) 

U(k = 2) 


0.010 
0.020 




0.60 

1.2 



Table 18c. Combined Standard Uncertainty for Thermal Resistance (R) forZ = 152.4 mm 













|c,- u\ 


(c,- u)ly 


*i 


x i 


c i 


u(x,) 




(m 2 -K-W _1 ) 


(%) 


Q 


0.8761 W 


-3.7599 m 2 -K-W~ 2 


0.0088 W 




0.03295 


1.0 


A 


0.12989 m 2 


25.3613 KW^ 1 


2.47xl0" 5 m 2 




0.00063 


0.02 


AT 


22.22 K 


0.14825 m 2 -W ' 


0.086 K 




0.01278 


0.39 


y = R-- 


= 3.29 m 2 K-W ' 






u c (R) 

U(k = 2) 


0.036 
0.071 


1.1 

2.2 



Table 18d. Combined Standard Uncertainty for Thermal Resistance (R) for L = 228.6 mm 













\Cj u\ 


(c t u)ly 


Xi 


x i 


C; 


u(x,) 




(m 2 -K-W"') 


(%) 


Q 


0.6112 W 


-7.72508 m 2 -K-W 2 


0.0084 W 




0.06454 


1.4 


A 


0.12989 m 2 


36.3531 K-W~' 


2.47 x 10 5 


2 
m 


0.00090 


0.02 


AT 


22.22 K 


0.21250 m 2 -W"' 


0.086 K 




0.01832 


0.39 



y = R = A.12m K-W" 



u c (R) 0.0068 

U(k = 2) 0.135 



1.4 

2.8 
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Using the values summarized in Table 18, Fig. 12 plots the relative combined standard uncertainty for A and indi- 
vidual uncertainty components for Q, A, and AT as a function of L. The graphical analysis is very similar to Fig. 10 
except for the absence of any uncertainty contribution from L in Fig. 12. The analysis again clearly shows that the 
contribution from Q is the major component of uncertainty, especially at thicknesses of 75 mm and above. The major 
component is the standard uncertainty of AT" which is fixed at 0.39 % across all levels of L. The uncertainty 
contribution due to A is nearly zero for all thicknesses. 
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Fig. 12. Combined standard uncertainty and individual components for R 
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11. Reporting Uncertainty 

The measurement result issued to a customer for a 
CTS specimen is for thermal resistance (R), not for 
thermal conductivity (A exp ). Therefore, the expanded 
uncertainty U and relative expanded uncertainty U r for 
thermal resistance (R) are reported as follows: 



Thermal Resistance (R): R, m 2 - K- W ' : 
m 2 - K- W \U n %) 



U, 



where the reported uncertainty is on expanded 
uncertainty as defined in the International 
Vocabulary of Basic and General terms in metrolo- 
gy, 2nd ed., ISO 1993 calculated using a coverage 
factor (k) of 2. The results given in this report apply 
only to the specimen tested, and not to any other 
specimen from the same or from a different lot of 
material. 

The relative expanded uncertainty U r for R is pro- 
vided to a customer with two significant digits that are 
rounded up to the nearest 0.5 %. For example, the 
values of U r for fibrous-glass blanket NIST CTS given 
in Table 18 are rounded as shown in Table 19. Values 



of U for the customer are rounded to be consistent 
with U r 

Table 19. Typical Values of R and U r for Low-Density Fibrous- 
Glass Blanket NIST CTS 



L 


R 


U 


(mm) 


(m 2 - K W" 1 ) 


(%) 


25.4 


0.564 


1.0 


76.2 


1.61 


1.5 


152.4 


3.31 


2.5 


28.6 


4.75 


3.0 



It is important to emphasize that other low-density 
fibrous-glass insulation materials may (and probably 
will) have different values of R and, consequently, 
different values of U r than those given in Table 19. 
Furthermore, it usually inappropriate to include in the 
uncertainty for a NIST result any component that arises 
from a NIST assessment of how the result might be 
employed [1]. These uncertainties may include, for 
example, effects arising from transportation of the meas- 
urement artifact to the customer laboratory including 
mechanical damage; passage of time; and differences 
between the environmental conditions at the customer 
laboratory and at NIST [1]. 
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12. Discussion 

NIST issues low-density fibrous-glass blanket CTS 
taken from an internal lot of insulation at nominal thick- 
nesses of 25 mm, 75 mm, 1 50 mm, or 225 mm. In gen- 
eral, measurements for customers are usually conducted 
at a mean temperature of 297 K and a temperature differ- 
ence of either 22.2 K or 27.8 K. Figure 1 3 plots the nom- 
inal thermal resistance (R) of the reference material as a 
function of specimen thickness (L) at T m of 291 K and AT" 
of 22.2 K. The error bars shown in Fig. 13 are equal to 
± U(k = 2) given in Table 1 8. 



Figure 1 4 examines the expanded uncertainties shown 
in Fig. 13 in greater detail by plotting U r from Table 18 
as a function of R. It is interesting to note that the trend 
for the expanded uncertainty data is non-linear. There are 
two possible related explanations for this non-linear 
behavior. At high levels of R, the specimen heat flow is 
reduced considerably thereby increasing the sensitivity 
coefficients (c,) and u c (Q) as shown in Eq. (32). Also, 
high levels of R are generally due to thick insulation 
specimens which will increase the edge heat flow effects. 



rr 



5.0- 




Error bars (± U) 




♦ 






P 


(kg-rrf 


3 ) 




4.0- 




9.4 
□ 9.0 










a 9.6 








3.0- 






8.7 




A 










2.0- 




□ 






1.0- 










n n - 


e 









200 



250 



50 100 150 

L, mm 
Fig. 13. Thermal resistance measurements of Fibrous-Glass Blanket NIST CTS (T m of 297 K, A^of 22.2 K) 
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Fig. 14. Relative expanded uncertainties (k = 2) versus thermal resistance of Fibrous-Glass Blanket NIST CTS. 
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Dominant Uncertainty Components 

Table 20 summarizes the individual contributions in 
percent (%) at the k = 1 level for A, L, AT, and Q 
(presented in Tables 17 and 18) for specimen thickness- 
es of 25.4 mm, 76.2 mm, 152.4 mm, and 228.6 mm. At 
25.4 mm, the dominant uncertainty component is for 
AT and, at thicknesses greater than 25.4 mm, the 
dominant component is for Q which increases consid- 
erably with L. The contribution of u (AT) is fixed for the 
measurements presented herein at 22.2 K; however, the 
contribution of u(AT) would change for different 
values of AT". It is interesting to note that the contribu- 
tion of L on A exp is the largest at 25.4 mm, decreasing at 
higher specimen thicknesses. Conversely, one would 
expect that, based on the results in Table 20, u cr (L) 
would increase for specimen thicknesses less than 
25.4 mm. In that case, u c (X,, xp ) could be larger than 
u c (R) because, strictly speaking, u c (L) does not directly 
affect the uncertainty computation for R as shown in 
Eq. (11). Finally, the small contribution for A, which is 
fixed for the all thicknesses, could change considerably 
at temperatures further from ambient due to thermal 
expansion effects (as noted above). 

Dominant Contributory Sources 

To investigate the results presented in Table 20 fur- 
ther, it is useful to re-examine the contributory sources 
for each component uncertainty of A, L, AT, and Q. 
Figure 1 3 reproduces Fig. 7 with the dominant contrib- 
utory source for each component identified (by circle). 
It is important to note that these dominant contributory 
sources may be different for other measurements 



(i.e., different specimen materials, operation tempera- 
tures, etc.). Furthermore, operators having other 
guarded-hot-plates will note that the dominant contrib- 
utory sources given here do not necessarily apply to 
their own measurement process. 

The contributory sources identified in Fig. 13 are 
discussed further below. 

• Meter area (A): At T m of 297 K and AT of 
22.2 K, u cr (A) is fixed at an estimate of 0.02 % 
for all measurements presented in this report. 
The dominant contributory sources are the 
uncertainties in plate dimensions of the guard 
gap (Table 3). However, at temperatures depart- 
ing from ambient, the contributory uncertainties 
for thermal expansion are expected to increase 
appreciably. 

• Specimen thickness (L): The dominant contribu- 
tory source for L is due to the uncertainty from 
the cold plate deflection under mechanical load 
(Table 6). This source is significant for com- 
pressible materials, such as thermal insulation 
blankets, and is typically smaller for more rigid 
and semi-rigid materials. In the case of rigid and 
semi-rigid materials, other contributory sources 
in Table 6 could become more important. 

• Temperature difference (AT): The dominant 
contributory source for AT" is the digital multi- 
meter (DMM) measurement of the PRT temper- 
ature sensor which has a standard uncertainty of 
0.058 K (Table 7). The standard uncertainty is 
based on the manufacturer specification and is 
probably a conservative estimate. 



Table 20. Percent Contribution of Individual Components for A and R 



25.4 mm 



A cxp (Table 17) 
76.2 mm 152.4 mm 



228.6 mm 



25.4 mm 



R (Table 18) 
76.2 mm 152.4 mm 



228.6 mm 



A 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


L 


0.15 


0.05 


0.02 


0.02 


... 


... 


... 


... 


AT 


0.39 


0.39 


0.39 


0.39 


0.39 


0.39 


0.39 


0.39 


Q 


0.17 


0.46 


1.0 


1.4 


0.17 


0.46 


1.0 


1.4 
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• Heat flow (Q): The dominant contributory source 
of Q is the uncertainty in the parasitic heat flows 
(AQ) across the guard gap and auxiliary insula- 
tion. The contributory sources for these parasitic 
heat flows are the control variability of the gap 
thermopile voltage of 1.5 uV or about 0.01 K (at 
the 3 standard deviation level) and the measure- 
ment of the PRT temperature sensors. 

Comparison With Previous Error Analysis 

Officially, the uncertainty assessment presented 
herein supersedes the previous error analysis prepared 
in 1983 [4]. This means that only the uncertainty values 
presented in this report should be used for the NIST 
Calibrated Transfer Specimens. Nevertheless, there is 
technical merit in discussing both analyses. Some of 
the obvious differences are the changes in measurement 
processes from 1983 to 2008 that involve different 
operators as well as modifications to some of the equip- 
ment (including both the apparatus and the data acqui- 
sition system). There is also a fundamental difference 
in the two approaches taken for the determination of 
combined standard uncertainty as discussed below. An 
examination of both approaches also reveals similari- 
ties in how individual and contributory uncertainties 
were determined. A brief discussion of the differences 
and similarities of the two analyses is presented below. 



Differences in approach: First and foremost, the error 
analysis prepared in 1983 [4] (hereafter, 1983 error anal- 
ysis) preceded the official NIST policy [1] that adopted 
current international guidelines for the expression of 
measurement uncertainty [2]. Consequently, the termi- 
nology given in the 1983 error analysis, as well as the 
approach for combining the uncertainties, has been ren- 
dered out of date. Other specific differences in the two 
approaches are given below. The 1983 error analysis: 

1) did not categorize the uncertainty components as 
either random or systematic as was done in an 
uncertainty analysis of an earlier NBS guarded- 
hot-plate apparatus [16]; and, 

2) estimated an "upper bound" of the "total uncer- 
tainty" by direct summation of the uncertainty 
components. 

In contrast, the GUM approach calculates a com- 
bined standard uncertainty by the root-sum-of-squares 
approach shown in Eq. (6). Further, the GUM and 
NIST policy [1] require that the expanded uncertainty 
(U) be reported with a coverage factor (k) equal to 2 for 
international comparisons. The 1983 error analysis 
does not report a coverage factor (k) and the "direct 
summation approach" makes it difficult to assess the 
coverage factor. Without a coverage factor (or a method 
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to determine one) it is difficult, if not impossible, to 
compare the combined standard uncertainties from the 
two analyses. 

Similarities in individual uncertainty components: 
A brief review of the individual uncertainties in the 
1983 error analysis revealed that several of the esti- 
mates agree reasonably well with the values presented 
in Table 17 for the combined standard uncertainty of 
A exp . Table 21 summarizes a side-by-side comparison of 
the individual uncertainties determined for this assess- 
ment and for the 1983 error analysis. 

For A and L, the standard uncertainties determined in 
this assessment and estimates of uncertainties for the 
1983 error analysis are in very good agreement. For AT", 
the standard uncertainties for this assessment are more 
conservative by a factor of approximately 2. As 
revealed in Table 21, the major differences in the indi- 
vidual components are the estimates for the specimen 
heat flow Q. The difference is attributed to the 
approach taken in this report for empirically determin- 
ing the uncertainty in parasitic heat flows of the meter 
area. The resulting Type B standard uncertainty 
u B (AQ), as observed in Table 16, dwarfs the other con- 
tributions for heat flow (Q) by, in some cases, two 
orders of magnitude. The author believes that the 
approach presented here for determining estimates of 
u c (AQ) using an experimental imbalance study is 
necessary to determine the standard uncertainty for the 
specimen heat flow (Q). 



13. Summary 

An assessment of uncertainties for the National 
Institute of Standards and Technology (NIST) 
1016 mm Guarded-Hot-Plate apparatus is described in 
this report. The uncertainties have been reported in a 
format consistent with current NIST policy on the 
expression of measurement uncertainty, which is based 



on recommendations by the Comite International des 
Poids et Mesures (CIPM) given in the Guide to the 
Expression of Uncertainty in Measurement. Strictly 
speaking, these uncertainties have been developed for a 
particular lot of low-density fibrous-glass blanket 
thermal insulation having a nominal bulk density of 
9.6 kg • m" 3 . This reference material, known as a NIST 
Calibrated Transfer Specimen (CTS), is issued to cus- 
tomers as specimens 610 mm square and at nominal 
thicknesses of 25.4 mm, 76.2 mm, 152.4 mm, and now 
228.6 mm for use in conjunction with the "representa- 
tive thickness" provision of the U.S. Federal Trade 
Commission (FTC) R-value Rule. The relative 
expanded uncertainty at a coverage factor of k equal to 
2 for the thermal resistance of this material increases 
from 1 % for a thickness of 25.4 mm to 3 % for a thick- 
ness of 228.6 mm. The approach for the assessment of 
uncertainties that has been developed herein is applica- 
ble to other insulation materials measured at 297 K. 

Recommendations for Future Research 

The uncertainty analyses given in this report have 
identified dominant components of uncertainty and, 
thus, potential areas for future measurement improve- 
ment. For the NIST 1016 mm Guarded-Hot-Plate 
apparatus considerable improvement, especially at 
higher levels of R, may be realized by developing 
better control strategies for the guard gap that include 
better measurement techniques for the gap voltage and 
PRT temperature sensors. In some cases, determining 
the individual standard uncertainties has required 
establishing traceability to NIST metrology laborato- 
ries, specifically for thermometry and primary electri- 
cal measurements. Recent calibrations from these 
metrology laboratories have indicated that more fre- 
quent checks and/or calibrations are required in the 
future. An extensive list of recommendations and future 
activities is given below: 



Table 21. Comparison of Individual Components in (%) for X„ 









A:xp 


(Table 


17) 






A cxp (from Table 1 


in Ref. [4]) 




X, 


25.4 mm 


76.2 mm 


152.4 mm 


228.6 mm 


25 mm 




75 mm 


150 mm 


300 mm 


A 


0.02 


0.02 




0.02 


0.02 


0.01 




0.01 


0.01 


0.01 


L 


0.15 


0.05 




0.02 


0.02 


0.1 




0.03- 


0.02 


0.01 


AT 


0.39 


0.39 




0.39 


0.39 


0.16 




0.16 


0.16 


0.16 


Q 


0.17 


0.46 




1.0 


1.4 


0.04* 




0.06* 


0.08* 


0.61* 



* Value obtained by summing heat flow estimates from Table 1 in Ref. [4]. 
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Annual check of platinum resistance thermo- 
meters (PRTs): The last calibration by the NIST 
Thermometry Group revealed a stability prob- 
lem with one of the cold-plate PRTs. To track 
any drift or potential problem, an annual check 
of the cold-plate PRTs at the triple point of water 
is now planned. 

Annual calibration of standard resistor: The last 
calibration by the NIST Quantum Electrical 
Metrology Division revealed a possible drift in 
the standard resistor used for the meter-plate 
power measurement. The presence of the drift 
could be indicative of detrimental factors affect- 
ing the resistor itself. Annual calibrations are 
now planned to investigate the extent of the 
possible drift. 

Improvement in thickness measurement: The 
dominant component for the thickness uncertain- 
ty is the analysis of the plate deflection under 
mechanical load. Because the approach present- 
ed here has limitations, an alternative technique 
to assess the plate deflection more accurately 
could reduce the thickness uncertainty. This 
would be useful at specimen thicknesses of 
25.4 mm or less. 

Improvement in temperature measurement: The 
dominant contributory source for the specimen 
AT is the digital multimeter (DMM) measure- 
ment of the PRT temperature sensor, which has a 
standard uncertainty of 0.058 K. The standard 
uncertainty for the DMM measurement is based 
on the manufacturer specification and is proba- 
bly conservative. Further, the uncertainty proba- 
bly includes systematic effects of unknown mag- 
nitude which will largely cancel when the speci- 
men difference is computed in Eq. (16). A signif- 
icant improvement in the temperature measure- 
ment could be realized by development of new 
instrumentation for measuring the PRT tempera- 
ture sensors. This improvement would also result 
in a lower uncertainty for the parasitic heat flows 
as discussed below. 

Reduction in the uncertainty of parasitic heat 
flows: The uncertainties of the parasitic heat 
flows are due primarily to the uncertainty in tem- 
perature measurement (and control) of the guard 
gap thermopile and the temperature difference of 
the auxiliary plate temperatures. As discussed 



above, reduction in the temperature measure- 
ment uncertainty for PRTs would significantly 
reduce the uncertainty in the parasitic heat flow, 
which is the dominant component of the speci- 
men heat-flow uncertainty. 

• Analysis of error due to edge heat transfer: 
Further work is recommended to investigate the 
differences between edge-heat flow errors com- 
puted from theory and those computed from an 
empirically-derived model presented in this 
report. 
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where: 



A mp = surface area of the meter plate (m 2 ); and, 

A gap = surface area of guard gap between meter 

plate and the guard plate (Figs. 1 and 2). 

For the circular meter area of the NIST guarded-hot- 
plate, Eq. (A-l) is rewritten as 

A = nr 2 +(nr -nr 2 )l2=-(r 2 + r) (A-2) 

" ' ° 2 



where: 

r = outer radius of meter plate (m); and, 
r, = inner radius of guard plate (m). 

Including thermal expansion effects, Eq. (A-2) is 
rewritten in the following form: 



A = -\r 2 (\ + aAT ) 2 +r 2 (\ + aAT )T 

2 L ° '"" ' mp J (A-3) 



^-[(r 2 +r 2 )(l + aAT mp ) 2 ] 



where: 



a = coefficient of thermal expansion of aluminum 
(alloy 6061-T6) (KT 1 ); and, 
A7L, = temperature difference of the meter plate from 
ambient (K) = r A -20°C. 

Further simplification of Eq. (A-3) can be obtained 
as shown in Eq. (A-4): 



A = 7tr 2 (l + aAT) 2 



(A-4) 



by substitution of: 



Appendix A — Meter Area (A) Derivation 

Test Method C 177 [7] defines the meter area as: 



A - A i s "> 



(A-l) 
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